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EXACT  SOLUTIONS  OV  S;iinT,E  VRT.OOrTY  KTHhTIC  fOPaWWt  AHD 


I 


THEIR  APPLIGATIOIv^  IN  CALCULf.TING  DIFFUSION  PROBLEMS 
(IMPROVED  DIFFUSION  METHOD) 
by  Yu.  4.  Romanov 

Kinetic  Equation.  Solution  of  T.'llnaea^s  Problem  and 

Determination  of  Albedo  of  Semi-infinite  Medium 

The  kinetic  equation  for  the  single-velocity  nlane  problem 
for  an  isotropic  scattering  Indicatrix  is,  as  is  Fell  knovm, 

in  rhich  n)  is  the  distribution  function  oi  the  number  of 
paT'ticles; 


ft  is  the  cosine  of  the  angle  betvr-een  the  direction 
of*  flight  of  a  p  .rticle  .nd  the  Z  axis; 

p  is  the  parameter  describing  the  properties  of  a 
medium  an'^  is  enual  po  the  ratio  of  the  scattering  cross  section 
and  the  total  cross-section  (p'<;.l). 


T’he  total  path  length  of  the  p.c.rticle  is  taken  as  unity. 


In  an  riedium  Eq.(l)has  the  exact  solr-tion 

'  y 

“  J  t(».  “  c, e**  -f  C,e-»*, 


-  I  - 


(?) 
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in  ’"'hich  the  particle  density; 


C  and  C  are  arbitrary  constants; 

“2 

k  is  the  coefficient  vdilch  is  found  from  the  solution 


of  the  transcendental  equation 


Afthil 

A 


t. 


(3) 


In  a  breeding  medium,  at  p  >  1,  the  roots  of  Eq.(_3)are  purely 
imaginary,  v-hUe  the  solutions  of  Eq.(^2)are  periodic. 

The  solutions  of  Eq. (2)satisry  the  diffusion  equation 


(4) 


provided  the  diffusion  factor  in  it  is  assumed  equal  to 


(5) 

(both  the  velocity  of  the  particle  and  its  path  length  are  equal 
to  unity).  Polutlon  (2)  describes  the  behavior  of  the  p'lrticle 
density  some  v.ay  from  the  interface  and  from;  the  sources  vith 
*.*■  f  ir  accur  ic'".  Eror  nov  on  the  oirticle  density  expressed  by 

Fq . '"ill  he  Icnown  -^9  *hG  asjrr.ptotlc  density,  v.-hile  the 

FTD-TT-61- 124/1+2 


2 


diffusion- coefficient  ^^is  knov.n  as  thc^ymptotic  diffusion 

coefficient.  Table  I  gJves  the  coefficients  k  and  for 

/O^ 

different  parameters  p  over  the  range  0.25  <  p  ^  2.0  . 

The  solution  of  Eq.(l)in  a  semi-infinite  medium  (iJilno^s 

problem)  can  be  found  by  the  Wiener  and  Hopf  method  /l/  and  is 

at*. 

given  in  /2~L/  for  ♦  particular  (of  a  non-absorptive  medium 


(p=  1).  Ref.  /5/  gives  a  solution  of  the  Milr^  problem  at  p  <C  1 
by  the  variation  method.  The  particle  density  function 
^d'-lno^s  nroblem  can  be  represented  as 


The  asymptotic  part  of  the  solution  l|^(from  novi'  on  called  the 
asymptotic  density),  as  can  be  seen  from  the  numerical  data  in 
/5/,  differs  from  the  true  density  t«(d)'  on  the  boundary  (z=0) 
by  loproxlmately  and  at  the  distance  of  the  semipath  from 


Lhe  bounlary  (7.=n.5'''  by  le.ss  than  3*^.  The  asyroptotic  density 

nCs)  in  Viln^^s  problem  is  e-jual  to  P*8e) 

*(*)  »!»*(*  +  *•). 


For  2  is  found  with  good  accuracy  from  the  following 

formula  ”  _  3  ^  1  ^  1 
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Table  1 


Values  of  coefficients  k, 


0 

k 

tK 

1 

/ 

k 

D* 

7 

0.95 

0.75102 

rp 

0.68 

0.84707 

0,44600 

0,7070 

**2S 

1*  ^  1 

0.74136 

0.69 

0,83804 

0.44148 

0.7176. 

0.t7 

1 '  1 

0.73179 

0.70 

0.82862 

0.43701 

0.7281 

0.98 

0.72231 

0.1651 

0.71 

0.81884 

0,43360 

0.7364 

0.18 

0.71294 

0.1810 

0.72 

0.80865 

0,42828 

0,7486 

0.73 

0.79804 

0.42402 

0.7£87 

0,00 

0.99741 

0.70364 

0.1969 

0.74 

0.4  IMI 

0,7187 

0.31 

0.99678 

0.09447 

0.2127 

0,75 

0.77551 

0.41M 

0.7786 

0,32 

0.08540 

0.2285 

0.76 

0.76354 

0.41157 

0,7885 

0,S 

H 

'  v;-?i 

0.07647 

0.3442 

0,77 

0,7982 

0.34 

0.06767 

0.3599 

0,78 

rG^Si 

UMJ 

0.8079 

0.35 

0,05898 

0.27.85 

0;79 

0,72454 

0.39998 

0.8174 

0.30 

1 SK 

0,65043 

0,2910 

0.37 

0,6420i 

0.3C64 

0.80 

0.39628 

0.8268 

0.38 

0.81375 

0.:'2I7 

O.Si 

0.8362 

0.39 

0.96743 

0.62563 

0.3369 

0.82 

0.68024 

0.84S5 

0.83 

0.66411 

0,8547 

0.40 

0.98*103 

4).6I763 

0.3520 

0,84 

0.84724 

0,38189 

0.8638 

0.41 

0,98364 

0.60979 

0.3669 

0.85 

0..17842 

0.8728 

0.42 

0,98150 

0,60307 

0.3816 

0.86 

0.61087 

0.37S0I 

0.8818 

0.43 

0.97917 

0.59451 

0.3962 

0.87 

0.59127 

0.37166 

0,8907 

0.44 

0.97667 

0.58707 

0,4I»7 

0.88 

0.57059 

0.3G839 

0.899SL 

0.48 

0.97397 

0.57979 

0.4251 

0.89 

0.548t8 

0.9083 

.0.40 

0.97109 

0,57363 

0.4.'<94 

0.47 

0,58562 

0,485 

0.90 

0..52543 

0,36306 

0.9170 

0.48 

I'K'!'.rrl 

0.55872 

0.4673 

0.91 

0,50061 

0,35897 

0.9356 

■■IIH 

0,92 

0.47397 

0,35'93 

0.0341 

0,49 

f  TTnYa 

0.S5I98 

0.41^08 

0,9.1 

0.44524 

0.35293 

0.9126 

0.50 

|i 

0.54537 

0.49-tO 

0,94 

0.41394 

0.51 

I*  v>is4J 

0..'3S8S 

0.5071 

0.95 

0.37948 

0,9593 

0.52 

■'  rtl’rif 

0.81252 

0,5.01 

0.96 

0.34081 

0,9676 

0,53 

|i  (’f|/(ylj 

0.526.’8 

0.7329 

0.97 

0.29625 

0,9758 

0.52016 

0.,5455 

0.98 

0,33868 

0,9839 

0.55 

0.9aVi3 

0.5I4I6 

0,5579 

0.99 

0,17227 

0.33598 

0.9920 

0.56 

0,50638 

0,5702 

0,57 

O.S«2SI 

0.5S24 

l.CO 

llACiXiS 

0.33.t'3 

i.oroo 

0,58 

0,49686 

0,5945 

1.01 

0.17,183 

1,(080 

0,59 

1 ' 

0.49132 

0,6064 

1,02 

0.24486 

pi 

1,0160 

1.08 

DBWl 

1,0236 

0.00 

0.90733 

0.48588 

0,6180 

1.04 

0.35188 

0.32315 

1,0316 

0.61 

0.90087 

41.4^055 

0.6296 

1.05 

0.39497 

0.32067 

1,0.193 

0.02 

0.89411 

0.47534 

0.6411 

1.06 

0.43431 

0.31822 

1,0470 

0.03 

0.88708 

0.47019 

0.6524 

1.07 

0.47092 

k0546 

0.04 

0.87973 

0,46510 

0.6636 

1.06 

0.85 

0,87306 

0.46023 

0,6746 

1,09 

0,53802 

0,31105 

1,0696 

0.08 

0,88407 

0,48539 

0.6855 

0.07 

0.85878 

0,45063 

0.8963 

1. 10 

0.56926 

0.30872 

1,0773 
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ii.3i 

’ll.** 

ll.34 

rtl.3T 

?J,3S 

Il,«> 

•1,44 
:  J.45 
1,45 
4,47 

JiS 

V,80 

1. 91 

t,tt 

1.91 

I, 94 

J, 8t 


%■ 


k 

Dt 

7 

1' 

O.S9926 

0,30641 

1,0848 

1,56 

0.b2H22 

0,30113 

1,0322 

1,67 

0,65676 

0,30189 

1,0996 

1.58 

0.6-951 

0,29908 

1,1069 

1,59 

0.7I00S 

0,.97S0 

1.1141 

0,73597 

0.29535 

1,1213 

1,60 

0,76135 

0,20323 

1,1285 

1.61 

0,76616 

0,29114 

1,1336 

1,62 

0,81056 

0,28908 

1,1427 

1.63 

1.64 

0,83454 

0.28705 

1,1498 

1,65 

0,65810 

0,26505 

1.IS68 

1.66 

0,88133 

0,2-309 

1,1638 

1,67 

0,90123 

0,28116 

1,1706. 

1,63 

0,92661 

0.2  925 

1,1777 

1,69 

0,94913 

0.27738 

1,1846 

0,97116 

0,27555 

1.1914 

1,70 

0.99299 

0,27375 

1,1982 

1.71 

1.01455 

0,27197 

1,2049 

1.72 

1,03593 

0,27021 

1,2116 

1.73 

1.74 

1,05707 

0,26816 

1,2183 

1,75 

1,07605 

0,2-672 

1,2250 

1,76 

1,09664 

0,26501 

1,2316 

1,77 

I.II947 

0.2602 

1,2382 

1.78 

1,13994 

0.26165 

1,2448 

1.79 

1,16014 

0.26000 

1,2513 

I,18*M2 

0,25836 

1,2578 

1,80 

1,20016 

0.25675 

1,26 <2 

1,81 

1,22036 

0.25516 

1,2706 

1,82 

1,24014 

0,25358 

1,2770 

1.83 

1.84 

1.25962 

0,25202 

1,2834 

1.85 

l,2?J37 

0,25049 

1.2897 

1,86 

1,29682 

0.24S97 

1,2960 

1,87 

1,31617 

0,24747 

1.30  >3 

1,88 

1,33742 

0,24590 

1.30S5 

1,89 

1,35658 

0.24452 

1,3147 

1,37565 

0,24306 

l,.1209 

1,90 

1.394ti2 

0,24165 

1.3271 

1.91 

1.410i3 

0,24023 

1,3312 

1,92 

1,43236 

0,23883 

1,3393 

1.93 

1.94 

1,49110 

0,23745 

1,3454 

1.95 

1,46978 

0.7306 

1,3514 

4.96 

1,48938 

0.23473 

1,3574 

1.97 

1,80603 

0,23339 

1,3634 

1.96 

1,82630 

0.332IM 

1.3094 

1.09 

1,94301 

0.2J077 

1.3754 

■ 

,  # 

3.00 

k 

t 

1,56216 

0.22948 

1*.3813 

1. 5-046 

0,22320 

1,3872 

1,59869 

0.22693 

1,3931 

1,61687 

0,22568 

1,3989 

1,63501 

0,22444 

1,4047 

1,65309 

0,22122 

1,4105 

I,b7|l2 

0,22201 

1,4163 

1,68910 

0.27081 

1,4221 

1 .7070'! 

0,21963 

1,4279 

1 ,72493 

0.21846 

1,4336 

1,74^79 

0,21730 

1,4393 

1,76059 

0,21615 

1,4450 

1,77835 

0,21502 

1,4506 

1,79609 

0,21389 

1,4562 

1,81378 

0,21278 

1,4618 

1,83143 

0,21168 

1,4674 

1,84905 

0,21059 

1.4730 

1.86663 

0,20951 

l.47te 

1.88419 

0.20H44 

1,4840 

1,90170 

0,20738 

1,4895 

1,91920 

0,20633 

1,4950 

l,93t«« 

0,20530 

1,5005 

1,95407 

0.2<M-27 

1,5059 

1,97146 

0,20326 

1,511.1 

1,98882 

0.20225 

1,5167 

2,00617 

0.20126 

1,5221 

2,02348 

0,20027 

1,5275 

2,01076 

0,19929 

1,5329 

2,05302 

0,19833 

1,5382 

2,075^5 

0,19737 

1,5435 

2,09246 

0,19612 

1,5188 

2,10961 

0,19548 

1,6541 

2,12661 

0,19455 

1.5594 

2,14396 

0,19362 

1,5646 

2,16106 

0,19271 

1,5691 

2,17815 

0,19181 

1,5750 

2,19524 

0,19091 

1.5BQ2 

2,21228 

0.19(712 

1..5854 

2,22932 

0,18914 

l.S!4)6 

2,24633 

0.18827 

1,5957. 

2,26333 

0,18740 

1.6008  ^ 

2.28011 

0,18655 

1,6069  ^ 

2,29725 

0,18570 

1,6110 

2,31422 

0,18485 

1,6161 

2,33113 

0,18402 

1,6212 
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in  i-hich  z,.  is  the  ex.tr -anolated  length,  for  v'hich  an  e-«'pres5ion 
is  given  heJov.  The  particle  flnx  on.  the  boundary  in  this  case 
is 


hrt'' 


y  w  - 1  •‘♦(o.  = - 1*(0)  J  i*T  (i*)  ♦•(0). 

tp-JOl  -  /  )  -  fcj 


The  angular  distribution  function  on  the  boundary 


itisfies 


the  equation 


_ L 


^jr+i*  '  2»(i*)(i -**!**)  ’ 


(9) 


Eqr. .  (_8)  anci  [9)  ir.akc  it  possible  to  obtain  a.-yTOximate  expressions 
for  the  function  <p(n)  .  In  the  .-^ero  y^upoxi tion  the  angular 

distribution  is  approxim.ated  by  *  linear  rational  function 

G 

and  the  const.anoe  a  and  b  are  found  from  the  normalization  condi¬ 


tion  and  the  exact  value  of  the  flux  on  the  boundary 


**  (see  previous  page) 

At  £  k  Is  Imaginary  and  the  solution  becomes  periodic;  the  author  aaria 

made  a  detailed  derivation  of  the  subsequent  equations  in  1951. 

-  6  - 


In  a  particul'Ar  case,  when  p=l,  Ep.  (lo)  becomes  the  angular 


distribution  formula  put  forv/ard  by  Fermi 


.  f(p) 


*+•? 


The  following  approximation  can  be  obtained  by  iteration  from 


Fn.[9) 


toll  _ 

The  function  Vi(l»)  has  -m  exict  O ;  at  other 

values  yK  the  accuracy  in  determining  9(11).  by  the  above  me;thod.  is 
0.1'^  . 


The  value  i*{p)  is  found  from  the  formula 


<*<»  —  1/  . — ! —  . 

y 


(11) 


rhe  function  »o(p)  can  be  represented  with  a  good  degree  of 


•iccur-:cy  by  the  formula 


(1^^) 


The  rroMe'"  of  refl'  ction  of  particles  from  an  infinite 
sn  ee  is  solved  by  the  '"iener-Hopf  m  thod.  we  find 

*-he  flbe-’o  '.nd  .ngul^r  di.c --pi^-utien  of  the  n  :rticles  esc  iping 
frc"  thi'  ’"f'dium  ♦-  (p)  from  the  'iven  .ngu’’  r  'irtriwution  of 

-  7  - 


incident  p,^.rticles  n  (^)  .  The  integral  relationship  linking 


♦+(l»)  and  t-(l»)(0<it<l)  takes  the  form 

t-W-yfWO 


(13) 


in  vrhich  the  function  f(l*)  i  s  the  angular . distribution  in  kiln^>s 
problem  and  satisfies  Eq .  ^9^  Equation ^13)  has  been  derived  by 
Ambartsuinyan  /6,7/  and  Chandresekhar  /8/  using  other  rethods. 

Using  the  integral  relationship  ^13)  v^e  can  calculate  the  albedo 
for  different  c  uses  of  incident  angular  distribution;  (a)  for 
a  flux  from  a  un'form  rlane  isotroric  source  situsAed  on  the 
boun  T  iry  ^t+(l»)=  . 


(1/) 


and  (b)  for  isotropic  distribution  of  the  incident  flux 

A  4-  f 

•  '*  '  /  (It) 

-.nd  (c)  for  -angular  distribution  of  the  flux  impinging  normally 
io  ^  he  Furf  ifiP  (t-f  (|()™A(}t— 1)], 


-  8  - 


(16) 


.S63iTjtl6n;;',of  Kinetic  Equation?  for  Tv/o  iwwitiiir  SfimMnflm’. t e 


Media 


Let  us  use  Eq.(l3)to  derive  a  system  of  inte'^ral  equations 
for  the  angular  distribution  on  the  plane  Interface  between  two 
lie  i.ia.  Let  the  particle  Ilux  be  directed  from  medium  1  to  medium  2, 
The  ^ingul'^.r  distribution  of  these  particles  will  be  designated  as 
K  (m)  ,  and  the  angular  distribution  of  particles  leaving  medium  ? 
-nd  impinging  on  medium  1  ’.vili  be  designated  as^.  (ji)  .  Since  on 
•iccount  of  the  conditions  of  the  oroblem  the  p-article  density  in 

eJcj 

n:edlum  ?  tends  to  zero/7 the  conns-ction.  between  the  functions 
♦♦  (1*)  'nd  (n)is  determined  by  Eq»  ^3) 


»-W=^  f  Wd  -  . 


(!?') 


In  v;.,loh  the  subscript  ,?  r:eans  thnt  the  correspen-'ln-  vajpee  end 
f'lnctlon?  refer  to  sedluit  ?.  The  other  e-iuatlon  Is  sl.Tllar  In  forn;, 
tut  It  r-ust  be  t-Jcen  Into  aooount  thot  the  particle  iensitT  in- 
ore  .r,es  exponentially  in  .  eilue,  1  .t  ^ 

nee  F-,.(;q«  have  to  seor-.te  fror,  U  (|.)-he  ^-.rt  of  the  angular 
nstrlbutionC»,(M)respon.-ltle  for  the  exponential  increase  in  the 

n*:rticli=>  ens’' ty  t  ^  Thus, 

-  9  - 


(17) 


yfiWO 

X  f*  9, (p.) rfH.. 

j  l‘+N 

in  xshlch  C  is  the  normalizing  constant  deteri’ining  the  flux 
coming  from  medium  1. 

The  solution  of  the  system  of  two  equations  (1?^)  and  (17) 
takes  the  form 


♦♦  (rt  =  ^Pt 


Ti(l*) 

(t  — *il«)Tt(V) 

(1— *iV)Ti((i)  ’ 


(IB) 


v.'hich  can  he  seen  by  direct  substitution.  The  normalizing  constant 


B  is  determined,  for  example,  from  the  condition  that  the  'density 
on  the  media  hounl'iry  is  equal  to  unity 


Tn  the  given  case  B  ~  while  the  flux  on  the  boundary  is  ex¬ 
pressed  by 


/  =  J (p) - 1- . 

'’’he  '■ngular  distribution  on  the  boundary  and. 
ca.lcul  •- ted  by  different  methods  in  /9jl0/. 

’"e  shoul  ’■'olnt  out  that  in  a  number  of  cases 


(19) 

the  flux  are 


eouations  of 


% 
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type  (l?^)  and  (17)  can  be  solved  if  the  solution  is  sought  in 


the  form 


n(,»)  =  /?,(,.)4. /?,(,»)  JEiM  , 

ftW 


(20) 


in  vjhich  are  linear  rational  functions. 

I'he  coefficients  in  the  linear  rational  functions  are  found 
from  algebraic  relationships  derived  after  substitution  of  ex¬ 
pressions  (?0)  into  the  initial  equations. 

The  particle  density  function^o(«) in  the  two  media  problem, 

as  in- ^lilno’s  problem,  can  be  represented  as 

^(x)  =  Cl-  #-*•*  h  Ci+  #*-»  +  ft  (z)  -  n^  (z)  +  /,  (*)  (z  <  0); 

Cj-  +/t(f)  !•/»(*)  >0)* 

in  ’Ahichfl(z)  is  the  asymptotic  density'-  and  f(7.)  only  differs 

from  ^.ero  in  direct  proximity  to  the  boundary.  The  asymptotic 

coefficients  and  are  found  as  deductions  with 

7l^  72\  /l-F 

r'-'-spect  to  the  ooles  s  =  +  of  the  Laplacian  mode  ♦!,*(*)  of 
the  function  ,  ’vhich  in  turn  is  expressed  by  the  angular 


distribution  on  the  boundary  iKO,  |i)  usin^'  the  foJ  lowing  formulae 

I  I 


»  J'  l+M  J  I— M 

- -  , 


I  -  -^Arth* 


f j1±M4_  _  f jizis 

/  t+p*  J  I- 


(?1) 


M 


-  i;-^Arlh* 

*  -  11 


The  results  of  Integfratlng  can  easily  be  shown  by  means  of  the  function 


and 


t.(*,  s)  = - — L - , 


(22) 


iOT 


t_(*.  s) 


«(*.*) 


f  W(*.  t*)4*  (*Vf  (*.  _  **  —  <* 

j  l+H  )  l-l«  I_i.Arth8  ■ 


c,_  s  lftL±Mid*iLiL>  , 

The  coefficients  ure  equal  to  .  *i) 

Cu  rs  *i)t(*i.  kty 

24|t— (Sj,  A|)t(A|8  A|) 

_  *(*!•  *») 
t(«^  ’ 

C»+  =0. 


(?3) 


The  solution  of  the  ■.■'roblen  corresponding  to  the  flux  directed 
frorr  nedlum  2  to  nedium  1  is  obviously  derived  by  replacing  the 
subscrlDt  2  by  1,  and  vice-versa,  in  the  first  problem.  The  general 
solution  rep-^esenls  a  linear  combination  of  the  two  solutions, 

"n-^  in  the  general  solution  the  asymptotic  density  satisfies  the 

* 

follov'ing  boun-^ary  conditions  ;  1.  FnuHli+-y  of  the  logarithmic 

at 


'derivatives  of  the  •’.symntotlc  density  in  the  extracolated  points 


Jl1£iL  =  j£lfiL 
•|(*|)  *|W 


(■/,) 


Zaretskly  as  well  (1952)  after  he  had  read  the 
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2.  Discontinuity  of  the  asymptotic  density  at  the  extrapolated 
points 


Wlfa)  -*»(*«)  , 

t(ri)  T(ft>/  (25) 

The  functions 

(P6) 

are  given  in  T'^ble  1. 

The  position  of  the  extrapolated  points  and  ^  in  Eqs.^24) 
ind[^25j '-^ne  found  by  the  follov/ing  formulae 

»,=/(*„  *,)-/(*„  K), 

in  rhich  z  is  the  distance  from  the  extrapolated  point  to  the 
interface;  'mrik  if  z.  is  positive,  the  extrapolated  point  is  in 
medium  2  (its  parameters  have  a  subscript  2),  whereas  if  it  is 
negative  then  it  is  found  in  medium  1.  The  extrapolated  length 
s  (Ec.  11)  is  also  linked  to  the  function  f,  to  v/it;  m  k). 

0 

The  function  f(k,s)  from  Ea.  5  is  determined  from  Sine  of  the 
equiv  ilent  formulae 
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(28) 


An  approKiination  of  Eq.  10  can  be  used  as  0>)  P<^)  is 
the  reciprocal  function  of  function  k(p)  determined  from  the 


transcendental  En.|3j. 


Table  2  shows  the  values  of  z.  To  determine  their  sign 


we  can  conveniently  use  thei  rule  that  both  extrapolated  points 
are  always  in  the  >  ore  active  medium  (in  the  medium  for  v.^hlch 
p  is  greater),  and  the  extrapolated  point  of  the  less  active 
matter  is  further  from  the  boundary.  All  the  preceding  derivations 
are  applicable  to  a  breeding  m.edium  (p  >  1)  and  the  corresponding 
formulae  hold,  since  k  is  purely  imaginary  in  this  case. 

To  complete  the  picture  let  us  write  down  the  relationships 
linking  the  true  particle  density  and  the  particle  density  on  the 
boundary  between  t\TO  media  with  the  asymptotic  density  and  its 


derivative: 


**(0)  —  I  /  — *!!_ 


iHiH)  . 
Tift)  ’ 


m 


K(i-ft)(i~ft) 


/*»*-*.»  ftlifL)- 
»(ft— ft)  Tift) 


ftlift) 

Tift) 


(29) 


.golutlon  of  Kinetic  rauation  V'ith  External  Particle  '-’ource 


T^suations  of  type and  (l7)for  determining  the  angular  dis¬ 


tribution  on  the  boundary  between  two  serai-infinite  media  can  be 


Distance  (z)  betv/een  extrapolated  point  and  interface  of  media 
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the 


At  2  for  the  adjoining  tnedlus  cont.  Table 
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used  to  solve  the  problem  of  the  particle  (Jistrlbution  around 


an  isotropic  source  located  on  the  boundary  of  these  media.  In 
this  case  the  equations  take  the  form 

kw=  ■  (30) 

I'he  source  on  the  boundary  is  normalized  in  such  a  way  that 


the  total  flux  from  a  unit  of  surface  is  equal  to  unity.  Equations 
(30}  are  solved  by  means  of  Eq.  20,  and  the  functions  K  (h)  and 

take  the  following  form 

A  *ii*)»i(rt  I  Pt  . 

^  2(fi— V)T*(l*)  2(ft— ^1)1*  ' 

Just  as  before,  the  oartlcle  density  is  represented  as  the  sum 

of  the  exponential  asymptotic  density  and  the  addition 


t,(»)  =  Ci+  +/,  W  (*  <  0); 

t#(*)  =  Ci-  +/,(ar)  (2  >0) 

iimw  On  iccount  of  the  condition  of  the  problem^  phere  is  no 


riux  from  inifinlty,  the  particle  density  at  z  and  z  — ^  “• 


tends  exuonentially  to  zero,  i.e.  =  0.  Just  as  before, 

/2+^ 

!he  asymntogic  coefficients  are  found  from  Eos.  21  and  are  eoual  to 
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•2*  ^ 


in  which  z 
"1 


C,+  -  k,)  _ 


¥ir: 


m-*,) 


^(Pt—P\) 


Ti« 


-Mi; 


’c«-  «  it.)t-(*..  *»)  _ 

0^  /l)  '  (4i»  *a) 


=  1/^ 


■V 


■Ml. 


(?2) 


tti  •♦•  kiUPt—pii 

and  z  are  the  distances  frd>m  the  extrapolated  points 


to  the  interfaces  determined  by  27. 


In  the  special  case  of /homof^enous  medium 


4p,(*)  =  .  g~k  111  4.  fi^\ 


and 


(33) 


This  partial  solution  is  contained  in  Karshak^s  review  /ll/. 

The  function  f(z)  is  positive  everywhere,  as  a  logarithmic 

anom-ily  at  z  =  0  and  decreases  ranidly  as  z  increases.  The  integral 

a 

of  function  f(z')  can  be  found  by  using  the  law  of  consetvhtlon oof 

v 

the  total  number  of  particles  (l-rt  J  t(*)*  -  i.  from  which 


1~30.T* 


folution  of  Kinetic  Equation  for  a  Cohere 
Vhen  solving  spherical  problems  in  which  the  particle  path 
constant,  it  is  very  useful  indeed  to  use  the  th'^orem 


* 


lencth  is 
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K 

of  the  reduction  of  a  sphefical  problem  to  a  plane  one  .  Accord¬ 
ing  to  this  theorem,  which  can  easily  be  derived  if  we  write 
down  the  kinetic  equations  in  the  form  of  spherical 

problem  with  a  set  distribution  of  the  function  p(r)  is  equivalent 
to  a  plane  problem  in  which  p(x)  =  p(-x)  =  p(r).  To  obtain  the 
particle  density  in  a  spherical  case,  we  select  an  odd  solution 
of  the  plane  problem  ifoW  ,  and  the  particle  density  in  the 
spherical  problem  ♦o(/-)  is  expressed  in  terms  of  the  particle 
density  in  the  plane  problem  by  the  formula 


♦.w 


(35) 


Here,  of  course,  the  angular  distributions  do  not  change  into 


one  another. 


This  theorem  implies  a  number  of  practical  corollaries; 

1,  In  an  infinite  spherical  problem  the  asyrnptotic  solution 
takes  the  form 


-ar 


(?6) 


'^he  asymptotic 


pl  ine  problem, 


density  in  the  spherical  problem,  Just  as  in  the 
is  the  solution  of  the  diffusion  eoU'-tion 
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with  a  corrected  asymptotic  diffusion  coefficient. 


2.  VTien  deriving  boundary  conditions  in  an  equi-path  spherical 
problem,  hr  has  to  be  substituted  into  the  boundary  oonditions  of 

•i*  ^ 

the  plane  problem  instead  of  function  n.  The  spherical  problem 
v.rith  a  constant  total  path  length,  in  which 


P(0=/»i.  .r<R  I 
r>R  \ 


reduces  to  a  plane  three-layer  problem 


(37) 


x<-R  h  x>R\ 

-R<x<R.  (og) 

Here  the  boundary  conditions  imposed  on  the  asymptotic  density 
on  the  boundaries  +  R.  are  the  same  as  in  the  problem  for  two 
semi-infinite  media  (Fqs.  2 A  and  25)).  Vdien  satisfying  the  condi- 
tlons  of  this  approximation  the  true  particle 

density  on  the  boundary  is  related  to  the  asymptotic  density 
by  the  following  equation,  easily  obtained  from  (^: 9), 


t.(«) 


V 


«l)  +  *1) 

R 


(39) 


derivation  of  the  flux  on  the  boundary  in  more  complex,  since 
^^ho  iheo^ei  c:'  re'uction  of  the  spherical  problem  to  a  pl'-.ne  one 
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when  the  riath  length  is  constant  is  not  applicable  to  the  flux. 


The  particle  density  function  may  be  represented  as  the  sum  of 
the  asymptotic  density  n(z)  and  the  f(z).  In  the  plane 


problem  the  particle  flux  on  the  boundary  (z  =  O)  is  expressed 


in  terms  of  the  integral  of  function  f(z)  ^  the  following  way 


y(0)  =  ~  -  (1  - /»,)J  {z)dz. 

On  the  other  hand,  Eq.  (^29) gives  us 


(40) 


T(fi) 


The  relationship  for  the  flux,  similar  to  Eq,(^40^'y takes  the 


form 


m  =  -  IV.  ^  J  /,  (r)  r-ir. 


(40') 

If  the  conditions  of  the  approximation  are 

satisfied,  by  virtue  of  the  theorem  of  reduction  we  obtalnt  the 
rel --tionshlp  a  .  If  the  total  path  length  of  the  particles  is 


constant,  by  expressing  f,-(r)  in  terms  of  f,(z'i  in  Fa.  40*^  ^.nd 

-slV 

substituting  the  integrals  of  f  If  (z)dz  is  found  from  Eqs.  40 

- 

fft(z)xdz 

and  :  vhile  /  is  calculated  by  differentiation  at  point  s=  0 
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« 


of  the  iaplifccian  modes  of  the  particle  density  in  the  plane 
problem  {/X)J ,  we  find  an  expression  for  the  flux  in  the  spherical 
problem 


Vfj?)-  -  _21|  4. 


(41) 


The  functions  f(k,s)  are  determined  by  Eqs.  2S;  in^4l'')'these 
functions  are  taken  at  s  =  0,  The  subscripts  i  and  e  ri.ean  that 
the  corresponding  parameter  refers  to  the  internal  or  external 
medium.  Eq.  ^l)  obviously  still  holds  when  subscript  1  is  replaced 
by  subscript  2  by  virtue  of  the  boundary  conditions. 

Use  of  Exact  Solutions  to  Calculate  Diffusion  Problems 
(improved  Diffusion  Method) 

The  approxim-ite  method  of  calculating  multi-layer  plane 
and  spherical  problems  for  media  with  a  constant  total  particle 
path  length  can  be  reduced  to  thr  folxoving:  an  asymptotic  particle 
density  function,  satisfying  the  diffusit^n  equation  with  an  asmmptotic 
diffusion  coefficient,  is  introduced  in  each  layer,  and  boundary 

I 

conditions  obtained  from] exact  solution  of  the  two  media  problem 
are  imposed  upon  it  at  the  boundaries  at  the  extrapolated  points. 

The  boundary  conditions  in  the  plane  problem/  Ets  ,  (p/.'^and  (?  5'^/  can 
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be  reduced  to  the  equality  of  logarithmic  derivatives  and 


discontinuity  in  asymptotic  density  at  the  extrapolated  points 
near  the  boundary.  As  is  clear  from  a  survey  of  American  research 
on'^inetic  theory  of  neutrons  /12/,  a  similar  method  called  the 
boundary  point  method  has  been  vforked  out  and  is  being  used  in 
American  calculations. 

®i 

Ouestions  of  the  accuracy,  limits  *prti%#oalil, li  I  j  and  [^ssibillty 
of  using  the  method  to  solve  spherical  problems  in  media  with  a 
varying  total  particle  path  length,  in  cylindrical  problems  and 
problems  with  distrlouted  sources  of  particles  can  best  be  illus¬ 
trated  by  solving  specific  problems. 

Critical  Dimensions""''^  ^ of  a  Plate  and  Sphere 

If  the  origin  of  the  coordinates  is  selected  in  the  center 
of  a  plate,  the  asymptotic  particle  density  in  the  plate  is  ex¬ 
pressed  by  the  function 


The  asymptotic  density  at  a  pwint^  the  distance  of  v/hich  from  the 
boundary  is  the  distance  of  the  extrapolated  lengthy  is  eounl  to 


sero 

*  For  Purposes  of  convenience  k  will  mean  the  Imaginary  part  of  the  trans¬ 
cendental  equation  root  from  now  on^t 


The  critical  thickness  of  the  plate  is 


*  (-J-  +*•)  =  «  co«  -  0. 


(4.2) 


We  should  point  out  that  the  same  result  is  found  if  we  use 
condition^.4)  for  the  equality  of  logarithmic  derivatives  at  the 
extrapolated  point  as  the  boundary  conditions.  Since  there  is 
no  reflection  we  can  assume  p  =  0  beyond  the  plate.  Hence,  k  =  1 
and  the  logarithmic  derivative  of  the  asymptotic  function  in  the 
ficticious  medium  beyond  the  plate  is  already  equal  to  -1,  Accord¬ 
ing  to  condition (2/^,  the  logarithmic  derivative  of  the  asymptotic 
density  In  the  plate  is  equal  to  -1  at  the  extrapolated  point 

1 .  e .  5 


At  p^=  0  /  at  p 


f  (t*)  “ 

0  ^  /,  Eqs .  (27yand  (r''S)give  us 


U'y<) 


•'fter  -‘lich  it  is  clear  both  critical  conditions(^42)  and  (4;^)are 
equivilent.  In  order  to  detern.ine  the  degree  of  accuracy  of  the 
formula  for  critical  d1m.ensicns  of  the  plate,  Table  2  compares 
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caloulations  of  the  critical  thickness  by  the  variation  method. 
4  formula  for  the  critical  dimensions  of  an  active  sphere  can 
be  obtained  in  similar  fashion 


ft 


.  0.71 


(43) 


/provided  z  is  represented  by  Eq«  1,2/. 
0 


Table  3 

Values  of  Critical  Thicknesses  of  the  Plane  Layer  Without  Reflector 


*  $  • . 

A 

-Up-  -  1 

‘  4  '  • 

‘  •  Mp  .  : 

rruui  uq* 

Variation 

method  | 

From  eq, 

W 

Varlat Ion 
method 

‘  0.1 

0.142 

0.157 

0,6 

’  0,482 

0.48 

0.2 

0.242 

0.2S1 

0,8 

0,560 

0.58 

0;4 

0,388 

■  0,39 

1,0 

.0,500.. 

.  0.50 

Comparison  ■'■"ith  nnnm'rical  calcul'dtions  shows  that  vhen 
determining  the  critical  radius,  the  inaccuracy  of  Eq.(^43)does 
not  exceed  1/  anywhere  over  the  range  of  p.  The  high  degree  of 
in  ccur  icy  in  the  given  case  is  due  to  the  ^sct  th''.t  the  critical 
r-.dius  -'-s  calcul'’t-ed  by^/3'jat  p  — >  co  chan-'-es  successfully  to 
limiting  forrul'i  ''•uoted  in  /I3/. 

Critical  hirrensions  of  tctive  Sphere  ’^urroundfed  by  Spherical 
L-ver  of  Inert  Reflector 

T^t  us  consider  -n  active  rph''re  '^'urroundtd  by  a  spherical 
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layer  of  inert  reflector.  Let  u.-  assume  that  the  particle  path 
length  in  the  reflector  is  eaual  to  the  total  length  of  the  par¬ 
ticles  in  the  sphere. 

The  asymptotic  density  in  the  active  sphere  (p  o  1)  is 
emjal  to 


Jli  •>».A 


mkr  • 

t 


At  r  =  0  the  function  n is  finite.  In  the  reflector  (r  =  1) 

The  following  boundary  conditions  are  imposed  upon  the  functions 

’-)  the  function  n,  disappears  on  the  external  extrapolated 
boundary  R  +  0.71 

7^ 


b)  on  the  boundary  of  the  active  sphere  (r  =  R)  at  the  ex- 
tr  isolated  points  the  logarithmic  derivatives  of  the  function  nr 


re  e-'uul  to 


4r  ,  ,  .  ■  ^  *r 

.V  t. ■ 


Ir  +  /( 
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Since  both  extrapolated  points  lie  inside  the  active  sphere 


(in  the  sphere  p  is  greater  than  in  the  reflector) , 


z  and  z 

^1 


are  negative. 


As  a  result  vie  obtain  a  formula  for  the  critical  dimensions 


=  71 (44) 

To  Illustrate  the  accuracy  and  limits  of  Eq*  44 


let  us  comp'^re  the  critical  radii  calculated  by  ^,.4)  and  calculated 


by  equation  for  p  =  1.724  (ic  =  I.B56,  z  =  0,0458,  z  = 

_  _  -  2 

-O.O865'',  at  different  values  of  the  relative  thickness  of  the 
reflecting  shell  *  r  ^  Table  4)  • 


Table  4 


Critical  radius  of  active  sphere  as  function  of  reflecting  shell 

thickness  (2.=  1.  724) 


‘S  ' 

r."'"  - — —  ■  — - — 

S 

From  equation 

from  numerical 

from  eq.uation 

ji  1 

calculations 

««0 

1.907 

1.981. 

s 

1.097  ' 

1.108 

— 

I.040 

1.0W 

0.99S 

1.008 

» 

0.955 

0.95« 

J 

0.937 

0.9*  ' 

0 

0.921 

0.98 

— 
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The  asymptotic  density  in  a  sphere  can  be  represented  as  the 
snm  of  the  homogeifeus  and  inhomogeneous  sdiutions  of  ikm  Eq.  30. 

On  function  n(r)  we  imposed  a  boundary  condition:/*^ /? 0  from 


which 


^  V  r,hik(/?+^)  ■ 


The  flux  through  the  external  surface  r  =  R  is  found  from  Eq.(^41^ 
As  a  result  we  obtained  the  multiplication  factor  0  v.'hich  is  equal 

tu 

to  ^  ratio  of  the  total  particle  flux  err.erging  through  the  external 
surface  and  the  total  oarticle  flux  of  the  source 


-P)  0) 

f  Sh*(^+*b) 


Ecu  ',tion(^45)  gives  us  the  mean  oath  length  of  the  particles 


in  the  sphere 


'=i‘.'!?XT  =  -«'  +  5^  +  T«+»W- 

'\t  P  1  ihe  error  in  Fq,(^/i6^doe?.  not  exceed  l'^-  ■=nd  at  P  =  0 

it  giver  an  incorrect  Umitin-  vaiue,  and  '.t  low  ve  can  obtain 


1  by  'be  '’uccessive  collision  rethod. 


^pnllcatlon  of  Method  To  Solving  Spherical  Problems 

IH  ^'edla^^hich  Total  Particle  Path  Length  Is  Not  Constant 
The  main  difficulty  in  solving  the  spherical  problem  for 
media  with  different  total  free  path  lengths  is  determining  the 
conditions  v»hich  hnve  to  be  imposed  on  the  asymptotic  density  at 
■^he  interface  .  It  has  not  been  possible  to  find  analyticalj^exact 
solutions  for  this  case,  but  the  derivation  of  boundary  conditions 
.  from  numerical  solutions  reauires  calculatlonlj^oundary 

conlditions  obtained  by  numerical  calculations  are  only  found  in 
scientific  litei-sture  for  a  few  individual  cases  (absolutely 
absorptive  spheres  and  cylinders  in  an  inert  medium) . 

In  one  of  his  research  projects  the  author  attempted  to 
compile  boundary  conditions  for  the  spherical  uroblem  with  a 
non'-constant  oartlcle  path  length  by  introducing  a  further  dis¬ 
continuity  in  the  neutron  density  on  the  interface  proportional 


to  the  flux  and  different  in  path  length.  But  the  boundary  condi- 

u 


tions  v.'hen  corrected  in  this  way  aimm-  a  small  sphere  of  applicability 


and  do  not  solve  the  set  task. 


The  fo]  lowing  as’^roxi'^ate  system  may  be  put  forv/ard  for 


calculating  the  eeneral  case  of  a  v.rying  path  length.  Fhen  there 
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are  no  sources,  the  follovdng  diffusion  equation  is  satisfied 
throughout 


div  / grad  +  (p  —  \)tn  =  0, 


(47) 


)h 


in  which  1  is  the  total  path  length  (or  transport-length,  AMTing 
non-isotropic  scattering,  in  the  laboratory  coordinate  syst'em). 
The  equation  contains  boundary  condition  for  n  on  the  disconti¬ 


nuity  boundaries  1  and  p: 


a)  continuity  of  the  flux  on  the  boundary 


/,  grad  -Om  «,  »/,grad 


(4S) 


b)  discontinuity  in  particle  density: 


i!^|  A  — 


(48^) 


Con''^itions  A  8^  ^  at  ^  =1  differ  froni  the  exact  ones  for  this 

1  ? 


c  s  e 


|lfc ondi t i ons ^ <  4^  ■'ind  ^ In  Eos .  ^4  ^Vand  (48^^  the  boundary  condi- 
ti  ons  (^Ipiposed  "t  the  true  bound  ries,  and  not  at  the  extrapolated 
points,  and  the  density  discontinuity  at  the  boundary  is  55il  , 


Am 


and  not  ^  .  Since  z  is  usually  considerably  smaller  than  the 

length  ,  and  the  function  is  almost  constant  ov-r  a  wide 
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range  of  variation  in  p,  this  approximation  proves  satisfactory 
in  a  number  of  cases.  Since  z  is  maximum  on  the  external  boundary, 
the  exact  boundary  conditions  are  best  kept  .  Instead  of 

conditions  the  follov-'ing  conditions  are  often  imposed  at  the 

boun^  ary 

% 

Ti  It 


These  conditions  are  not  a  corollary  of  the  differential  equation 
?.iut  are  more  exact  in  a  number  of  cases.  It  should  be  pointed 
out  that  in  h  case  in  W'hich  the  path  length  in  the  outside  layer 
of  the  spheric. -.1  .system  tends  to  infinity,  condi tions  (^/j 8)and  (>,.83 
■do  not  give  a  correct  limiting  process. 

Problem  Kith  Distributed  Source  in  Sphere 

To  solve  the  problem  using  theWthoery  of  pertubatlons  ,  let  us 
replace  the  true  distribution  of  sources  q(r^  by  an  eouiv.alent 
(f)  -  Htir)  the  value  ^  being  found  from  condition 

tyti'dV-  (;,q) 


.fter  this  the  problem 


fin  tjww  the  critical  value  of  the 


P  .rame^er  ’//  -p+5  at  a  known  value  of  the  sphere  radius  R  /from 

the  eigenfunction  in  the  critical  stite.  From 

*  In  the  general  caae  of  varable  Pelerls  parameters  and  distribution  of  the 
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» 


isotropic  sources  ^  ,  the  averaging  conditions  for  the  theory  of 
per'turbation  take  the  form 

.  .  .. 


(f  is  the  density  of  the  particle  flux).  The  condition  for  averaging  the 
parameter  is  written  down  in  the  diffusion  approximation.  The  theory  of 
perturbation  for  a  kinetic  equation  was  developed  by  Fuchs  and,  independently, 
in  the  USSR  by  N,A.  Dmltriyev  (1948), 
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the  condition  of  particle  balance  v/e  find  the  outward  flux  and 
the  multiplication  factor  which  is  equal  to 


Q  =  =  I  - 

(50) 

N  state 

is  exact  near  the  critical  *xtK,  but  its  accuracy  is 

satisfactory  in  other  cases  as  well,  even  at  p  =  0  (absolutely 

absorptive  inedium) . 

For  a  source  distributed  throughout  the  eigenfunction,  Ea.  (^5o") 
is  exact  for  any  values  of  R,  and  the  multiplication  factor  is 


Q=l  + 


Below  we  ?ive  formulae  determining  the  eigenfunc-tions  and 
its  integrals  for  ».  case  of  an  active  sphere.  The  integral 
in  the  critical  state  is  expressed,  according  to  the  perturbation 
theory  ,  in  terms  of  the  deriv^itive  of  the  critical  radius  with 
resu'ct  to  the  parameter  p 

The  ’tnee'^r'-:!  ^  is  fouiri  the  ’■'■'rtiele  nurber  pil'ince;  thus, 

•*  ,  .  -  — 

d  /-I 
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For  an  active  sphere  in  the  critical 
the  follov.'ing  relationships 


give  us 


il5L=:v^£±Z<*L0L , 

^  *  ft 

l-i'or  a  source  vi’ith  a  constant  density,  we  obtain  from  p]n.(50) 


Q=l  + 


3Do(p')lp-l) 


^  /p 

For  a  central  point  source,  Eq.(^5o) gives  us  a  forn’ula  coinciding 
near  the  critical  state  wdth 

This  method  if  calculation  v/as  then  applied  to  the  solution 
of  a  number  of  specific  problems. 
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A,p-prny1giite  Method  of  Calculating  Critical  Masses  of 


_  Spherical  Reactors  with  Infinite  Reflectors 

by  G.  I.  Marchuk  and  V.  P.  Kochergln 

Let  us  consider  a  critical  reactor  in  vliich  the  active 

■.one  of  the  r-'dlu.s  R,  imii  which  is  a  spherically  symmetrical 

neutron  source,  is  surrounded  by  an  infinite  reflector:  there 

the 

are  no  neutron  sources  in/reflector.  If  the  origin  of  the  co¬ 
ordinates  is  selected  in  the  center  of  the  active  zone,  the 
solution  of  the  single-group  diffusion  eauation  for  the  neutron 
fluT  in  the  reflector  takes  the  rorm  /!/ 


The  effect  of  the  reflector 

of  the  reactor  is  described 


'♦(r)  =  C 


(1) 


on  re-iucing  the  critical  dimension, s 

OK 

by  the  effective  additiw 


In  which  extrapolated  radius  of  the  reaccor  vithout  .a 

r- fleeter . 

"he  e:'f»^ctive  id-liticn  can  -also  be  defined  as  the  length  of 
^he  line- r  ex'r  relation  of  the  nantron  flux  in  the  reflector  on  the 
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boundary  v.'lth  the  reactori’s  active  zone 


iR 


- =-JL 


2. 

R 


fror/i  the  equality  of  Eqs .  (2)  and (3)  v/e  can  determine 


(3) 

ar 


X- 


R,^R  R 


(4) 


If  and  R  y'^y^are  knovai,then  by  solving  Fn.^4y».ith  respect 


to  R  re  obtain 


2X 


(5) 


To  fn.-..ke  Lhe  rndlns  positive,  -i  plus  sign  is  put  before  the  root. 

I  convenient  "fuantity  describing  the  energy  spectrum  of  the 

wan 


reflectorless  reactor  is  the  so-called  cadrnii 


chamber  ratio 


•J-’w 

CdR -  jyt/f  (u)itu/  jyl/f  (u)clu.  (6) 

••tu  -•> 

V  m 

The  value  A  of  1  he  ‘’unction  of  the  cadr.iihM:--. tio  is  determined 

•y  .  '■■•ns  Urn  of  the  critical  ogrameters  and  R  ^  derived  by 

solving  the  mult-igrour  dl ''fusion  enuntlons  at  the  'iivcn  densities 

of  ihc-  nncle  -r  fuel  nd  rolerator,  as  veil  as  at  the  given  uranium 

'■35 

Gnrich'"cn1  v-llh  resrect  +0  isotope  U  (cor  ex'u'.plo,  lOOb) . 
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I  < 

Fig*  1.  Function 

UOi-HiO. 


for  spherical  reactors. 


The  critical  size  of  .i  reactor  with  an  infinite  reflector 


with  other  enrichments  and  the  same  fuel  and  moderator  densities 
is  calculated  from  Eq«^5)nsing  ^^graph  for  the  function  \  .  But 


first  a  '.’lultigrouo  calculation  of  the  reflector 


reactors  has 


to  he  made  to  determine  cadmium  ratio. 

\ 

“^he  <tv  ^vh3  af  the  I'unction  K  for  •^yrterns  v;ith  vrater,  i^T‘;irhlte 
anda**  reflectors  are  riven  in  Figs.  1,  ?  and  Irrespec¬ 

tively.  Comrarison  of  ‘he  critical  dlrensions  of  the  -.ctlve  zone 
for  th  se  rv~  i  e^-s^  d' terrrlned  by  r'e  >ns  of  the  multigroup  calculation 


In  < ’ip  ■  rpT’s  ay  -..rschuk  .n'’  others  (see  nn.  ;-9,  and  57  in  the 
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!l!!!lll!millBBIIII  ■iimiii 


4  ttMfi,  *  •  ttn-'t  *  $ $1,0 1  4  $$»  /  « 

Cadmium  ratio 

»  ■ 

. _ 

Pig,  2,  Function  R^  —  R  R  for  spherical  reactors  U-C, 
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Critical  masses  of  uranium-plutonium-water  systems 
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✓ 


-41  . 


the  same  as  in  Table 


present  collection')  an' 
in  the  ratios«  =  p^Pu«, 


-|l5v'  ^  YTL-iA_gJL 

by  E<?.  (^5)iwease  vitl.e  ayhep*  of  variation 


are  viven  in  Tables  1  and  ?.  The  asterisk 


denotes  the  radius  found  by  (5). 

Discrepancies  in  the  critical  rass  lie  within  the  limits  of 
accuracy  vrith  which  the  multigroup  calculation  w^as  made,  and 


do  not  in  effect  exceed  10)o. 

Since  \  is  determined  by  means  of  the  multigroup  calculation, 
it  thereby  takes  into  account  the  retardation  effect  of  the  neutrons 
in  the  feflector;  hence  X  can  be  both  positive  and  negative, 
according  to  the  effectiveness  of  the  reflector. 

It  follows  from  Eq.  5  that  vhen  V  varies  bet\reen  -  cao  and  +  oo, 


R  v'ries  between  0  and  R^,  and 


The  value  \  retnrins  p 
6  g/cm  )  im  the  of 

consi  ' er  hion,  -*ni  the  eff 
th*n  H /v/b  throughout;  At. 


ositive  for  the  system  UO*-  H,.0  = 

A 

variation  in  the  cadmium  ratio  under 
ectlve  addition  is  therefore  small 

a  certain  value  ^he  cadr-iur-  r-’tio- 
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er 


•ff 

o- 


I 
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in  the  systems  U  -  C  and  U  -  Be,  \  changes  Its  sign  and  becomes 
^legative,  im±  by  virtue  of  which  the  effective  addition  in  this 
region  becomes  greater  than  R  /^,  This  also  follows  from  the 


Tables . 


iTo 


Te  should  point  out  that  inf^lculation  of  the  critical 


dimensions  of  the  systems  (df  •=  T  =6  g/cm  )  from 

dd  1)  VvOi. 


En.^5')it  is  nossible  to  use  the  function  X  shown  in  Fig.  1, 


confirmed  by  Table  1,  in  which  K  ,  and  the  uranium  is  natural. 
This  fact  shov.s  that  the  graph  of  oho  function  X  does  not 

change  when  uranium  is  replaced  by  plutanium^  provided  the  compound 

in  which  they  are  present  exhibits  the  same  density.  It  probably 

re¬ 
does  not  change  eithei"  when  uranium  or  plutonium  are/placed  -by 

.nothcr  fissll^  element.  Thus,  ’"e  .are  in  a  position  to  calcul  te 

the  critical  masses  :^rom  ET.^5)for  any  combinations  of  fissile 

elements  1«xrt)h  in  the  glv°n  comnoahd,  if  the  graph  of  function  X  is 

■=vail  .Me . 

~o  :eter:  in--'  the  critical  ’i.'.sses  (>f  s<‘'hfcrical  reactors^^  uiiK 
‘nfini^e  reflectors'  t  o'^her  fuel  ind  moder  .  cor  iensiti'^s,  it  is 


■c 


essenti  il  to  Ci7.culate  the  critical  naram''-ters  't'-y  means  of  inulti- 


:^rcur  '.Pions  Mr  ii  ff  usior^ 
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any  enrichment  of  the  fuel 


<VT=W  / 

fto  determine  the  function  V.  . 


For'''^irRXKtiar  reflectorless  reactors  there  exists  a 


/3/,  tfeK  and  the  critical  dimensions  can  easily  he  converted  when 
“he  fuel  and  moderator  densities  vary.  Let  us  now  consider  a  system 


conslstinf;  of  two  components. 


If  v;e  fix  ®  "Pf/Pjn.»  the  number  of  fuel  nucle»  per  1  cm  is 


determined  in  the  followin'-^  w-ay 


fc- 


-T  +«"4- 


(8) 


in  which  p*«  _^y  (au  is  the  Avogadro  number; 
g/cm  ;  A  is  the  xtornic  v.elght) . 


is  the  density, 


If  and  'I'  vary,  Pf  and  Pm  w'ill  vary  in  identical 

proportion,  hence  to  ch  -nge  from  one  system  to  the  other  the 

- ~  5  =  fp/j5[f» 

rnicroscoeic  auiirg  sections  have  be  multiplied  by'  JL  *'  .  Here, 

so  as  not  to  -'isrurt  the  multigroup  'Iffuslon  equations,  the  linear 
^i"ensions  h 've  to  be  divided  by  the  sii  llarity  factor  ^  . 


'i’hus  , 


^=-r 


(9) 
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The  curves  sho^rlng  the  critical  mode  as  a  function  of  the  ctitical 


volume  when  and  change  in  the  same  ratio  are  transposed 

■f  "  fn 

parallelly  v'ithout  changing  the  shape,  since  the  similarity  factor 
^  in  this  case  is  the  same  for  all  (t  . 


If 


and -V  do  not  change  in  the  same  ratio,  and  the  similarity 

**'  i  ^  o'Cu-tosQ. 

factor  therefore  depends  on  (X  ,  the  luBinl  Tninmr  curve.^are  still 


luiiil  liiiirim  the  transfer. 
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Footnote  on  page of  English  text: 

In  physical  meaning  K  is  MhMMc  the  reciprocal  of  an  effective  diffusion  length. 
Some  authors  use  \  in  the  effective  bovindary  conditions  method  /2/ .  In  this 
Yfork  \  is  considered  from  a  different  vie-wpoint. 
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OF  EVEN  '.PPPOXIU.^^TIONS  IN  SPHERIC 4L  HAKJONIC  I..ETHOD 


by  G.  Ya.  Rurnyant'^ev 


^f^hen  solving  the  kinetic  Boltzmann  equation  by  the  spherical 
harmonic  method  v^e  make  use  of  the  so-called  approximations, 
in  v.hich  N  is  the  order  of  approximation.  In  particular,  the 


vfell-knorn  diffusion  theory  is  identical  to  the  approximation. 

'a 

If  it  is  necessary  to  make  the  diffusion1Si»  theory  more  accurate, 


v-'o  can  resort  to  Hporoximations  of  higher  orders,  and  it  is  usually 

understood  tint  this  means  the  R.  and  P.  approximations,  or,  in 

"A  ' 

4 

other  vords^  approxim  ations  :«c  exclusively  uneven  orders.  Ro  far 
there  h-’.s  been  no  mention  in  scientific  l.lter-'ture  of  any  exa-imples 
of  the  pr  ctical  .ipnlication  of  the  an eroximations  of  ven  orders. 

Tht-  possibility  of  using,'^^  auproxiiaatlon  is  of  particular 
Interest  since  there  is  reason  to  believe  ::h'it  It  ought  to  be  more 
accurate  than  the  P.  •■ipuroximation  fsince  it  is  one  o."  the  I’ollovln'^ 
order 'I  ;ne  .t  the  s  ime  time  less  cumnersome  than  she  apprdxiivi- 


Men.  In  "'--ny  :--ses  it  could  he  s--.tl''f-ctorily  used  xz  x  cor-i.’sctL 


on 


to  the  ;if ‘uslon  theory.  Instead  of  P  . 


"he  reison  vrhy  the  even  •»epvoxim  ttons  -n’e  ieno'^ed  Is  th-.t 
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the  nunber  of  geniiine  solutions  in  them,  and  conseouently  the 
number  of  arbitrary  coefficients  in  the  .general  solution,  is 
smaller  than  the  number  of  spherical  components  retained  in  the 
e^cpanslon  of  a  distribution  function.  On  account  of  this. it  is 

, _ _ _ _ t, 

Impossible  to  satisfy  [on  the  interface  of  tv;o  unlike  media/ those 
boundary  conditions  which  are  usually  imposed  in  the  case  of 
uneven  aporoximations,  to  vit,  conditions  of  continuity  of  all 
these  spherical  hammonlcs  individually,  (it  should  be  pointed 

out  that  in  nont-unidimensiono.!  problems  this  difficulty  also  arises 
in  the  c:..se  of  uneven  ap-'roxim.’-tions . )  Belov  ve  L’ive  a  brief  descrip 
tlon  of  this  thod  of  formulating  the  boundary  conditions,  enabling 
us  to  use  the  "'.D'oroximations  of  any  orders,  including  _P^  and  the 
other  even  ones.  For  the  Sake  of  simplicity  we  illusL.r..te  the 
method  with  the  Rr^mple  of  a  plane  unidimenslonat  nroblem.  '.s 
r-pg  ;ris  symbols,  we  are  keeping  to  D.  Galanin'’s  book  /l/. 

Let  us  condider  the  system  of  equations  for  the  spherical 
harmonic  method 


*  irri/.-i + 5^/;.! = 0  (0  < « <  „), 


(1) 


This  system  contains  an  Infinite  number  of  equations.  The 
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that  In  the  equation  corresponding  to  n=N^ 

Is  discarded,  anA  the  system  then  becomes 
solution  In  a  form  such  that 

(2) 

we  reduce  the  differential  equations  to  normal  algebraic  ones. 

From  the  condition  that  the  determinant  of  the  system  is  equal  to 
zero  (condition  of  compatibility),  we  find  for  a  characteristic 
equation^  the  roots  of  which  are  ^genvalues  of  the  problem. 

It  can  be  shown  that  the  characteristic  equations  for  different 
p^approximations  will  take  the  form  of  polynomials 


P  approximation  is 

tf+i 


the  term 


finite.  Seeking  the 


Ar=i  «*-3*=o. 
iV-2  (5  +  4t)«*-I5*=:0, 

AT.- 3  9a«-(55€  4.35)a»+ I05i»0, 

AT  =  4 .  (64e  4-  I6I)««  -  (735t  +  315)a*  +  945.  =  0. 

Here  *  —  1  —  »  =  -p  . 


O) 
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This  type  of  equation  corresponds  to  a  case  of  isotropic 

scattering,  when  all  the  c  ,  except  for  c.=  1,  are  equal  to  zero, 

'n  7b^ 

This,  however,  is  wm  of  no  importance  in  principle. 

As  we  see,  the  degree  of  the  characteristic  polynomials  with, 

respect  to  CX  is  equal  to  N  +  1,  if  N  is  odd,  and  N,  if  N  is  even 

(this  regularity  exists  as  well  for  N  ::=s  4)  •  Consequently,  the 

number  of  ^g^values  of  0^’^ »  and  therefore  the  number  of  partial 

solutions  determined  by  Eq.^2)at  different  oij  will  always  be  even  ^ 

and  equal  to  N  +  1  or  N,  according  to  the  parity  of  the  approximation. 

Thus,  the  number  of  arbitrary  coefficientts  in  the  general  solution 

only  coincides  with  the  number  of  spherical  hamroonlcs  for  tmeven  N, 

for  in  the  P  approximation  for  a  plane  unidimensional  problem  the 
'N 

number  of  them  is  equal  to  N  +  1. 

Eqs.^^l)  still  hold  even  when  their  coefficients  reflecting  the 
properties  of  the  medljm  are  arbitary  functions  of  the  coordinates. 

We  can  therefore  consider  that  they  have  been  written  down  for  the 
medium  as  a  whole,  amd  that  in  a  heterogeneous  medium  the  material 
characteristics  1,  ^  and  c  take  the  form  of  plecewlse-constant 
functions. 

Let  us  integrate  each  equation  from  the  arbitrary  lower  limit 
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> 


a,  and  get  as  a  result 


m 

■•■[irr,/— +  5^a/.*iJ,.,- 


(4) 


There  is  always  a  function  under  the  integral  sign,  hence 
the  right  hand  side  of  Eq.(4)is  always  continuous.  This  implies  the 
requirement  of  continuity  in  the  expressions 


(5) 


This  requirement  should  be  satisfied  as  well  on  the  Interface 


between  two  media. 


The  combinations  J  aril  coefficients  of  the  serial  expansion 

7n^ 

with  respect  to  the  spherical  harmonics  of  the  function /(x,  «)cos^. 


Indeed, 


aw-f  I 

4« 


J/co.*/>,(co.»)<«  = 


(6) 


Consequently,  the  system  of  conditions  (5)  can  be  considered  as 
the  continuity  condition  of  the  function /(x,<>)  cosO.  .  As  we  see, 
conjdltions  (5)  are  equivalent  to  the  continuity  requirement  of 


the  integrals 
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J/c08*P*(COS*)rfQ  (0</»<oo). 

At  mk  0  we  have  a  continuity  condition  for  the  dlffusloh  flow. 

In  the  P  approximation  a  ,  which  implies  the  continuity 
“N 

•f  .  Taking  this  fact  into  account  and  using  J  ,  we  arrive  at  ^ 

V  -/  'n-1 

■econtinuity  f  .etc.  It  is  not  possible  to  draw  a  similar  conclu- 
~N-l' 

Sion  with  regard  to  the  functions  f  ,  f  ,and  so  on.  They  are 

~N  “N-l' 

only  part  of  the  continuity  condition  in  the  form  of  combinations 
(5). 


Thus,  conditions  on  the  boundary  between  the  two  media  reduce 
to  the [condl tion  of /^ontirmTEy jtir  the  following  values. 


/tt-i;  fn-i:  fit-%  •"  /wo; 

/v-i;  Jh-%  /hi). 

(8) 

The  minimum  value  of  the  subscript  in  functions  (^s)  depends  on  the 
parity  of  IT.  The  number  of  Joining  equations  is  obviously  even  in 
every  case.  It  is  ^difficult  to  calculate  that  for  an  uneven  N  this 
number  is  N  +  1  ,  and  for  an  even  N  it  is  N,  l.e.,  it  also 
coincides  with  the  number  of  arbitary  coefficients  in  the  solution. 
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In  particular,  for  the  approximation  we  only  find  two  condl- 

^2 


tions,  expressed  in  the  continuity  of  the  functions 


/i(jf)  »  /•(•«)+ y/jW- 

Typical  of  even  approximations  is  the  fact  that  the  scalar 

flow,  which  coincides  with  the  function  f  (x)  up  to  an  accuracy 

-0  - 

of  the*iiii  normalizing  multiplier,  experiences  a  value  at  the 

Interface  whftfth  drops  as  N  Increases.  When  there  is  no  absorption 

in  either  mediun,  the  discontinuity  disappears, 

in 

For  uneven  approximations  in  problems/whtch  the  number  of 
conditions  (s) coincides  with  the  number  of  spherical  harmonics,  the 
Itm  conditions  derived  are  Identical  to  Ik  the  continuity  condition 
for  each  separate  harmonic. 

It  has  been  demonstrated  in  Ref.  /2/  that  in  the  general  case 
conditions  of  type  ^8)should  be  derived  from  continuity  on  the 
boundary  of  the  InteKrals 


J/(ra)cos(a*)K^(2)rfa(-«<m<«,  0<s<AO. 
the 

in  which  V  is  normal  to  /vhboundary. 


(9) 
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Y^^are  spherical  junctions  of  the  general  type  which  can 
be  selected  ,  for  instance,  in  the  form 


K..*.=P*,(cos»)/  7*'”^ 
(  slnmtp 


0/ 

are  the  Leg^dre  joint  polynomials. 

(See  /3/  and  bther  courses  in  mathematical  analysis  for  spheri¬ 
cal  functions  in  greater  detail.) 

Conditions  (9),  the  ntimber  of  which  in  the  general  case  in  equal 
to  N(N  +  1),  make  the  problem  a  closed  one,  irrespective  of  the 
natxire  of  the  geometry  or  the  parity  of  the  approximation. 

The  numerical  examples  and  theoretical  arguments  show  that 
in  these  cases  in  which  the  use  of  approximations  of  a  low  order 

is  generally  possible,  the  P  approximation  can  make  the  P  approxi- 

"2  '1 

mation  considerably  more  accurate.  At  the  same  time,  the  laborious 

nature  of  the  P  approximation,  compared  with  the  P  ,  Increases 
'2  ^1 


only  slightly  in  most  cases  of  practical  interest 
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CRITICAL  MASSES  OF  UE.'UaiM  GRAPHITr  REACTORS 


BY  G.  I.  Marchuk,  G.  A.  Ilyasova,  V.  Ye.  Kolesov, 

V.  P.  Kochergin,  L.  I.  Kuznetsova  and  Ye.  I.  Pogudalina 

Introduction 

No  detailed  infonnation  is  available  so  far  in  scientific  literature  on  the 
critical  masses  of  uranium  graphite  reactors.  In  view  of  this  the  need  arose  to 
make  the  relevant  calculations  for  reactors  v/ith  a  wide  variety  of  neutron  energy 
spectra.  Some  very  interesting  calculations  have  b-’cn  made  by  Safonov,  hut  they 

O  "5 

relate  to  reactors  with  100:"  enrichnent  vdth  tho  uranium  isotope  /l/ . 

This  paper  deals  vdth  the  problem  of  the  critical  masres  of  uranium  graphite 
reacto-’'s  vdth  different  degrees  of  enrichment. 

Basic  Equations 

In  most  cases  the  diffusion  approyimation  is  suf^dcient  to  calculate  the 
critical  masses  of  reactors.  It  is  assumed  that  the  neutrons  in  tlue  reqetor  are 
inodrrated  by  elastic  and  non-elastic  scatterr  g.  I'ho  effect  of  neutron  bhermaliaation 
is  taken  into  account  in  the  low-neutron  ‘"•nergy  region. 

Pne  multi  group  system  of  equations  Tor  th*--  reactor  takes  the  form 
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E/^sr  — y/->  —  jt/Q(r) 

.  ^'  =  I/t' 

.1.1 

Q(7)^£i£}  vS/f  ®  +  Qif 

/•I  ' 


and^ 

moderation/'-'non-elastic  .gcattering,  and  the  total  removal  and  fission  cros,s-sections, 
while  D  is  the  diffusion  coefficient 


in  which  Sj  %  is  the  ela.stic  3cat^ering  cross  section,  and 
Nfc  is  the  mean  cosine  o'"  the  scatr.'rin:  angle. 
In  ’ir:-  hirh-energj'  region,  it  is  advisable  to  average  the 
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scattering;  capture  and  fission  cross-sections  with  consideration 
for  the  fission  spectrum,  and  in  the  intermediate  energy  region 

X 

with  consideration  for  the  Fermi  spectrum  ,  The  constants 
and  J2/»Tfor  the  heat  group  have  to  be  averaged  over  the 

steady-state  neutron  spectrum  in  an  infinite  homogeneous  medium^ 

« 

taking  into  account  the  thermal  motion  of  the  moderator  nucleie. 
For  this  purpose  we  use  the  model  of  a  single -atom  gas  moderator 
put  forward  by  Koen  /2/.  The  corresponding  equation  for  neutron 

XX 

density  N  derived  in  Wilkins*  differential  form  /2/  is 


(2) 


V7kT 


in  which  (v  is  the  neutron  velocity,  P  is  the  a1»solute  tem¬ 

perature  of  the  medium,  and  k  is  the  Boltzmann  constant); 

is  the  moderating  power  at  E  =  0.025  evj?. 

The  function  N(x)  at  x  =  0  satisfies  the  condltlonc'N(O)  ^  N*(0)=0^ 
ir'(oy=  const.  Vfhen  solving  EQ.(^2)we  find  the  function  N(x)»  which 


is  then  used  to  average  the  physical  constants: 


/f? 


N{x)ix 


^xN{x)ix 


V  r  *rp 


H{x)dx 


I 


xN(x)dx 


in  which  is  the  boundary  separating  the  thermal  diffusion 

region  from  the  moderation  region,  and 

£^,  is  the  absorption  and  fission  cross  section  for  thermal 

neutrons  /3/. 

Considerable  attention  should  be  given  to  the  resonance  neutron 

capture  in  the  intermediate  energy  region.  Twelve  resonance  levels 

235  ; 

have  been  calculated  on  me  U  nuclew  in  the  capture  cross  section 

23S  ' 

and  four  levels  on  the  U  nucle4n.  The  corresp6ndlng  probabilities 
of  avoiding  the  resonance  capture  for  a  particular  isotope  are 
calculated  by  the  formula 


<T>*  =  e  ^ 
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Here  n  is  the  resonance  number,  and 

fh  is  the  number  of  nuclei  of  the  absorber  per  unit  of  volume. 
The  effective  resonance  integral  is  determined  by  the  formula 


-aIL*L 
'V I  +  *« 


in  which  ir  is  the  total  resonance  integral  at  the  level  with  the 
number  n,  and 

T  is  the  temperature  of  the  medium; 

A  is  the  absorption  cross  section  at  maximum  re- 

•HI 

sonance; 


V  is  the  potential  scattering  cross  section  ); 

.i|-  is  factor  taking  into  account  the  Doppler  broadening 

of  the  resonance  lines  due  to  the  effect  of  the  thermal  motion  of 

r  ' 

moderator  nucleM^; 

.1,  —  — is  the  resonance  width; 

*  ' 


V 


The  function f|(t  A)  is  turbulated  in  A/. Markelov  and  Tyutere^ 
have  put  forward  a  simple  Interpolation  formula  for  calculating 


this  function 
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in  which 


ji  +  /F+ItTI  -VF-Pt+i.  «  o,09<«<i. 

l*.05.  '  if  {>1. 

p=A(0.15  +  0;  ' 

TP  =  l!p_  0,2285 +  0,282. 


This  formula  was  the  one  used  in  the  calculations. 

If  the  group  with  number  j  contains  so-called  resonances, 
we  calculate  the  total  probability  of  avoiding  resonance  capture 
in  the  group  <y>.  Here  the  density  of  the  moderation  of  the 
neutrons  iipe  in  the  group  J 

v=i:r;' 

l-l 

has  to  be  multiplied  by  the  factor  <ip^>  , 

It  should  also  be  noted  that  there  is  a  definite  probability 

235 

of  fission  at  the  resonances  of  U  .If  the  probability  of  resonance 

capture  with  fission  is  desifa*ted  1/(1+  then  the  probability 

of  radiation  captWee  is  equal  to  ^=^(1+  c^).  The  total  number  of 

235 

fissions  caused  by  resonance  capture  of  neutrons  on  U  nuvlel  takes 


the  form 


in  which 


m 


» -  <»>« 
!+•-  .  ’ 
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and  summation  is  with  respect  to  all  the  U  resonances  reaching 
the  group  with  number  j. 

Multigroup  System  of  Constants 

The  multigroup  of  physical  constants  has  been  compiled  on  the 
basis  of  published  data. 

the  fission  spectrum  region)  a  system  of  constants  has  been 

238 

J  irffief.  /5A 


235 

put  forward  for  U  and  U 


In  the  region  of  inter¬ 


mediate  energies  the  system  of  constants  was  derived  on  the  basis 
of  experiraental  daW  processed  by  Malyshev  /6/.  The  constants  for 


the  heat  group  are  taken  from  Ref.^^.  The  averaged  group  constants 
for  graphite  have  been  checked  by  calculating  the  neutron  age  in 


the  graphite.  The  ne»tron  age  up  to  the  indium  resonance  is  318 

2  -fe 

cm  ,  which  tallies  well  with  the  experimental  value .\The  multigroup 

system  of  physical  constants  is  given  in  Tables  l-3>  and  the 

235  238 

resonance  parameters  for  U  and  U  are  given  in  Tables  4  and  5, 
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TABLE  2 


Matrix  of  Inelastic  scattering  ^ 


/ 

* 

1 

2 

3 

4 

5 

6 

a 

o.'so 

3 

0.49 

0.38 

-  i 

4 

0.35 

0.45. 

.0.38 

ar 

s  • 

.0.14 

0.22 

0.10 

0.14 

.6 

0.05 

0.09 

0.U3 

0.03 

0.076 

7 

0.02 

0.03 

O.OJ 

0.01 

0.01 

0,025 

TABLE  3, 


Matrix  of  Inelastic  scattering 


y 

1  '  ■ 

1 

2 

. 

3 

4. 

5 

6 

1 

2 

0.71* 

• 

i3 

0,60 

0,40 

*  •• 

•• 

4 

0.50 

0.35 

0.34 

•• 

— 

0.20 

0.14 

0 

0.24 

4 

0.07 

0.05 

0 

0 

0,15 

0,02 

0,02 

0 

0 

0 

0.068 

Matrix  of  Inelastic  scatterldj^ 
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9'  O 
1 

« 

e 

•r 

/ 

E 

c 

1 

11 

35.3 

1730 

9.35 

0,987 

0.4 

2  ' 

II 

34.6 

1510 

5.76 

0.692 

0.69 

3 

11 

33.7 

800 

3,72 

0,835 

0.50 

4 

11 

32.2 

700 

3.40 

0.851 

0.50 

5 

11.25 

23.6 

900 

6.00 

0.996 

vO.50 

6 

11.50 

21.1 

-  700 

5.05 

1.06 

0.50 

7 

11.50 

19.40 

1590 

15.4 

1,32 

0,62 

8 

12,00 

12.40 

'2310 

18,4 

0.854 

1.7/ 

9 

12.00 

11,65 

1840 

lO.S 

0..*96 

3*33 

10 

12.25 

8.82 

1970 

33.8 

1,57 

0.67 

II 

12.75  : 

6.39 

1120 

14.5 

1.02 

0,96 

12 

13.00 

4,84 

530 

5.0 

0,630 

6,25  . 
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Table  5 


23s 

Resonance  Parameters  for  U 


The  multigroup  system  of  equations  for  reactors  without  re¬ 
flectors  takes  the  form 

!•! 

«=£'<V^T'+yt/.* 
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(3) 


.  (  /  <f 

17 

Here  in  groups  In  which  there  are  resonances,  and 


in  the  remaining  groups; 

is  the  geometrical  parameter  and  the  first  Eiganvalue 
of  the  problem 


U) 


is  the  extrapolated  reactor  surface) . 


It  is  assximed  that  In  the  resonance  energy  region  the  modera¬ 
tion  density  q  is  expressed  in  terms  of  the  flux  by  the  formula 


The  problem  is  solved  by  the  method  of  successive  approximations, 
which  is  normally  called  the  method  of  source  integration  /?/. 
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•  t  4  4  f  r  /  4  4  »«*  }  t  s  *4f*  t  **»«•;  *  #  •**  *  t  Hit 

Critical  volume  n 

Fig.  1.  Critical  masses  of  spherical  uranium- graphite 
reactors  without  reflectors  (y^  ■  1.65  g/cm^) 

Broken  lines  join  points  with  Identical  a; 

Solid  lines  show  different  degrees  of  enrichment 
with  Isotope  (%). 


We  should  point  out  that  transition  from  an 

extrapolated  surface  to  a  true  surface,  the  extrapolated  length 
was  averaged  over  the  energy  neutron  spectrum  in  the  reactor. 
SoMtlon  of  the  problem  gives  the  critical  loads  and  critical 
volumes  of  the  spherical  reactors  without  reflectors  mm  a  wide 


range  of  the  ratio  at  different  degrees  of  enrichment 


of  the  uranium  by  Isotope  U  . 


The  results  of  calculation  of  the  critical  masses  of  the 


reactors  are  given  in  Fig.  1 
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f  I  i  -ASin^t 

-.  Critical  volume, 4 

Fig.  2.  Critical  masses  of  uranium-graphite  reactors  without 
reflectors  as  function  of  volume. 

Solid  line  -  theoretical  data  ,  uranium  enrichment  90%. 

■  1.65  g/cm-^ 

o  -  experimental  data,  uranium  enrichment  93.2%  y  -  1,645  g/cm^ 

c 


Fig,  2  compares  the  results  of  the  calculations  and  measure- 
*.!ents  made  in  uranium-graphite  reactors  without  reflectors  /8/. 
The  close  correspondence  of  the  experimental  and  theoretical  data 
is  observed  in  the  region  of  both  thermal  neutron  reactors  and 
incermediate  neutron  reactors. 
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Fig.  3  shows  the  cadmium  ratio  for  the  fission  chamber  of 


I 


the  function  of  a.  ,  The  curve  shown  corresponds  to  a  lOOJ^  uranium 


enrichment.  The  curves  for  other  uranium  enrichments  colncidd/luor 
practical  purposes  with  the  ones  in  shown  Fig*  3. 


Figs.  4  5  show  graphs  for  minimum  critical  masses  and 


their  volumes.  Pig.  6  compares  the  results  of  calculation  of 
wrltlcal  masses  with  experimental  data  for  pure  uranium  as  a 
function  of  enrichment  /9/.  Analysis  that  the  diffusion  approxima¬ 
tion  leads  to  substantial  errors,  even  for  small  systems. 
Calculation  of  Spherical  Reactors  With  Reflectors 
Reactors  with  graphite  reflectors  used  to  be  calculated  by 

("A 

tie  multilgroup  method,  the  principle  of  which  is  described  in 
Ref,  /?/.  It  was  assumed  that  the  reflector  was  of  infinite 
thickness,  but  that  the  calculation  was  made  with  a  reflector 
thickness  of  70  cm.  The  error  due  to  substitution  of  a  finite 
for  an  infinite  reflector  was  comparable  with  the  errors  of  the 
numlerical  calculation 
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Fig.  7.  Economic  advantage  of  graphite  reflector  for 
spherical  uranium-graphite  reactors  Ye  »  1,65 
Figures  on  curves  show  enrichment  with  isotope 


(%) 


0u> 

Fig.  7  gives  the  effective  additions  aef  a  function  of 
at  different  degrees  of  enrichment^  derived  imtias-  numerical  solu¬ 
tion  of  multigroup  reactor  equations.  The  maximum  value  of  the 
effective  addition  for  reactors  with  a  100^  uranium  enrichment 
is  found  in  thermal  reactors  and  is  equal  to  the  square  root  of 
the  migration  area  iR  =  64  cm,  and  decreases  monotonically  to 
8  cm  for  fast  reactors.  For  intermediate  neutron  neactors  we 
observe  an  increase  in  the  effective  addition. 

Given  other  degrees  of  uranium  enrichment,  the  curve  showing 
the  effective  addition  of  a  function  of  oC  is  also  complex  in  form. 
When  the  enrichment  decreases,  the  maximum  effect  of  addition  is 
attained  in  the  region  of  Intermediate  neutron  s pec trum%ac tors. 
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It  may  be  assumed  that  there  is  a  region  in  which  infinitely 


large  suberitical  reactors  without  reflectors  may  become  critical 
if  reflectors  are  used.  This  means  that  an  infinitely  multiplying 
medium  coinciding  in  properties  with  the  active  zone  of  the  reactor 
is  suberitical  (k  ^  -<1).  At  the  same  time,  if  the  size  of 

the  active  zone  surrounding  the  reflector  is  restricted,  the  system 
may  become  supercritical. 


It  should  be  pointed  out  that  to  calculate  fast  neutron  reac¬ 


tors^  the  diffusion  approximation  is  not  sufficiently  accurate.  Fig, 
8  gives  the  critical  masses  of  spherical  uranium-graphite  reactors 


with  infinite  graphite  reflectors. 


Fig.  8.  The  crltlcfl  aasMt.of  uranium- graphite 

reactors  with  Iniintte  grapnite  rerlectdrs.  Broken 

llnaa  loin  ftpiata  with  ideatlcai  solid  lines  show  e'nrlchment 
with  Isotope  U***  (X), 
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CRITICAL  MS8ES  OF  URANIUM-BERYLLIUM  REACTORS 


by 


K 

G,  I.  Marchuk,  G.  A.  Ilyasova,  V.  Ye.  Jfolesov,  V,  P.  Kochergin, 
L.  I,  Kuznetsova  and  Ye. I.  Pogudalina. 


A  number  op  papers  /!/  deal  with  the  calculation  of  critical 
masses  of  uranium  -beryllium  reactors,  but  they  only  touch  on 

calculations  of  reactors  without  .refle6‘tors  with  lOOJ^  uranium 

235 


enrichment  with  U  .  The  present  paper  contains  the  results  of 

calculating  critical  masses  of  uranium-beryllium  reactors  both 

without  (reflector  and  with  an  infinite  beryllium  reflector  at 

235  4- 

different  degrees  of  enrichment  v-ith  Isotope  U  •  Alithe  cal¬ 
culations  have  been  made  in  the  diffusion-age  approximation,  by 
the  multigroup  method  Het  forth  in  the  paper  by  Marchuk  and  others 


entitled  "Critical  masses  of  Uranium-graphite  Reactors",  in  which 


uranium-graphite  systems  are  investigated  in  a  similar  fashion 


Group  constants 


Letargy 

•2? 

aB« 

• 

/-/ 

/Ib* 

~I.63S>- 0.375  . 

0.44 

1.5 

0.050 

0.35 

0.375-1.0 

0,65 

2.6 

O.O'iO 

0,08 

1.0  -1.6.'S 

1.00 

3.4 

1,000 

0.02 

1. has -2.25 

0.90 

3.9 

0.900 

2.25  -3,0 

1.10 

4.8 

1,100 

3,0  -3,75 

1.30 

5.4 

1.300 

3.75  -4,5 

1.3 

5.5 

1.300 

4,5  -5,25 

1.66 

5,60 

1.66 

5,25  -7.0 

0,710 

5,56 

0.710 

— 

7,0  -8,5 

0,828 

5.56 

0.828 

8.5  -9.5 

1,24 

5.56 

1.24 

— 

9,5  -10.5 

1.24 

5.56 

1,24 

— 

10,5  -11,5 

1,24 

5.56 

1.24 

— 

11,5  -12.5 

1.24 

5.66 

1.24 

12.6  -13.5 

1.24 

5.56 

1.24 

— 

13,5  -14,5 

1.24 

5,55 

— 

14.5  -  15,5 

1.20 

5,46 

1.20 

— 

Note;  for  the  heat  group  »J|*=rO,OI;  (I  —  co* »)“*  6,«;  t»5J  71.46. 


Fig.  1.  Critical  masses  of  spherical  uranium-beiTrllium  reactors  -.rithout  reflectors 


tot  4  ttto^  t  TJJWt  4  t»tO*Z  Vtt/Q* 

Critical  voliane,  1, 

Fig,  2,  Critical  masses  of  spherical  uranium-berj'-llium  reactors  vdthout  reflectors 


g/cm^),  i’he  meaning  of  the  lines  is  as  in  rig,  1, 


Enrichment,  7. 


Fig.  3»  Minimal  critical  masses  of  ui’anium  beryllium  reactors  ('Y 
1)  Tdthout  reflector;  2)  v/ith  beryl'' ium  reflector. 


1 . 8  g/ crn^ ) . 
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Fig.  U 


Fig. 


Critical  masses  of  uranium-beryllium  reactors  with  beryllium  reflector 
1,8^  g/cm^).  laeaning  of  lines  as  in  previous  figures. 


Zf’ective  additives  for  reactors  vri.thout  reflectors  as  ftnction  of 
for  three  enrichments  (V  -  1*5?  g/cm^). 
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In  the  high  energy  region,  apart  from  elastic  scattering. 


we  calculating  nonelastic  scattering  in  uranium  and  the  reaction 

(n,  2n)  on  beryllium.  The  corresponding  group  constants  are  given 

in  Refs.  /3/.  To  obtain  the  group  constants  in  the  intermediate 

energy  region  we  used  data  from  Ref*  /3/.  Here  the  calculations 

235  238, 

took  into  account  the  resonance  absop^tlon  on  U  and  D  xivt 
u  and  at  low  energies  took  into  account  the  thermal! zatlon  of 
nettrons  ^vlthin  the  framework  of  the  model  single-atom  gas  moderator 
/4/,  The  heat  constants  for  beryllium  are  taken  from  /5/. 

The  multigroup  system  of  constants  for  beryllium  is  given  in 
the  Table.  The  age  of  the  neutrons  up  to  the  indium  resonance  energy^ 

calcula-^lllpiiHM.  using  these  group  constants,  is  81  cm  (beryllium 

^  3 

dansity  Yg  =  1.85  g/cm  ),  and  this  tallies  well  with  the  experiment, 
The  results  of  the  calculations  of  uranium-beryllium  reactors 
without  reflectors  at  ratios  «=^PBe/pu“  (in  which  ^  is  the 

nuclear  d«i41ty)  and  fllfferent  degrees  of  enrichment  with  isotope 

235 


U  are  given  in  Figs,  1  to  3.  The  critical  load  curve  correspond¬ 
ing  to  a  100/b  enrichment  tallies  well,,  qualitatively  speaking,  with 
the  data  in  /!/. 

When  calculating  reactors  with  a  reflector,  the  thickness  of 


the  latter  was  taken  as  49  cm.  A  beryllium  reflector  of  this 

thickness  can  be  considered  infinite.  Results  of  the  calculations 

235 

of  the  critical  masses  of  U  and  the  effective  additions  are 


given  in  Figs.  4  and  5  respectively. 

Calculations  were  made  with  a  Strelai^ (arrow)  computer. 
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CRITICAL  MASSES  OF  WATER  MIXTURES  OF  URANIUM  AND  PLUTONIUM  COMPOUNDS 


by 

G.  I.  Marchuk,  G.  A,  Ilyasova,  V.  Ye.  Kolesov,  V.  P.  Kochergin, 
L.  I.  Kuznetsova  and  Ye.  I.  Pogudalina. 


Knowledge  of  the  critical  masses  of  different  wm  uranium  and 
plutonium  compounds  is  very  important  both  in  designing  and  building 


nuclear  reactors  as  well  as  in  solving  problems  safety. 


Systems  with  wmImmammImmiitKUat  a  water-containing  rr.oderator  are 
of  great  practical  Interest.  The  results  of  some  experimental 
investigations  of  critical  masses  of  uranivim-«ater  mlii%ures  have 
recently  been  published.  In  particular.  Ref.  /!/  discusses  a 
large  number  of  experiments  made  by  scientists. 

It  is  an  extremely  difficult  thlnjj  to  calculate  the  critical 
masses  of  uranl\iia-water  and  plutonium-water  reactors.  The  main 
difficulty  is  that  hydrogen  moderation  cannot  be  considered  con¬ 
tinuous,  and  on  that  account  the  age  theory  of  calculation  proves 
Inapplicable.  So  far  comparatively  little  theoretical  data  on  the 
critical  masses  of  water-containing  systems  has  been  published. 


We  can  therefore  only  refer  to  Sasonov's  calculations  which  are 
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given  in  /2/ 


The  authors  made  a  large  number  of  calculations  of  the  critical 
masses  of  reactors  toth  without  reflectors  and  with  water  reflectors. 

The  fuel  used  was  UO  with  water  and  a  mixture  of  UO  +  PuO  with 

2  2  2 

water.  The  concentration  of  the  fissile  matter  and  the  mixture 
varied  within  wide  limits.  The  calculations  were  made  for  different 
degrees  of  enrichment  and  different  ratios  of  the  number  of  uranium 
nuclei  to  plutonium.  The  error  in  the  critical  mass  did  not  exceed 


All  the  calculations  were  made  with  a  Strela  liomputer^ 
Basic  Equations 


The  multigroup  system  of  reactor  equations  in  the 


diffusion 


approximation  in  the  presence  of  hydrogen  can  be  as  follows 
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V  V 


vW  +  +  x'Q(r). 

+  Ti. 


$mt 


•  it 


I  ■  t—m+t 


(1) 


Qir) «  E  Vt  +  Qk  . 


/-I 

y  /i-y+i  i-y+i\ 

^ "  E  (  Si«  +  E*  j  To  "I"  •••  **) 

{U\_ 


t-j 

Here  Ii.  ‘ 


are  the  group]nonelastlc  scattering)  cross  sections^from 


the  group  with  ntunber  1  into  group  v/ith  the  number  Jj 

are  the  group  moderation  cross  sections  on  all 
nuclei,  except  hydrogen, 

fill  are  values  describing  the  elastic  scattering  on 
hydrogen  nuclei  from  group  1  into  group  Jj 

is  the  total  cross  section  of  from  group 

The  heat  constants  were  aver4|ged  by  the  foriaull 


tHx)dx 


1 


mIHic)Mjc 


k 
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~yN(xydx 

—  Rot)  + 


in  which  —  I^t)  are  the  thermal  cross  sections  at  an 

energy  0.025  ev, 

is  the  boundary  separating  the  thermal  diffusion  asfc 
region  from  the  moderation  region, 

T  is  the  temperature  of  the  medium, 
is  the  room  temperature. 

The  function  N(x)  is  the  neutron  spectrum  in  an  Infinite 
homogeneous  medium^allovdng  for  the  thermal  motion  of  the  moderator 
nuclei. 


To  find  N(x)  we  use 

y'(x)  -  S{x)y{x)  -  R(x)y{x) 


with  the  initial  conditions 

•  y(0)  =  0,  /(0)  =  const. 
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Here 


5(x)  — P(x)l/l:erf(x); 

R(x)  =  P(jf)  X*  +  r{x)  —i—  ; 

K«Pfx) 


The  resonances  in  the  absorption  cross  sections  were  taken 

into  account  in  the  following  way.  In  the  groups  with  resonances 

\ 

the  moderation  density  was  multiplied  by  the  probability  of 


avoiding  resonance  capture  <f'>  •  The  calculation  of  was 

done  in  the  same  way  as  described  ln|  ’'Critical  masses  of  uranium- 


graphite  reactors"  by  iiiarchuk  and  others  ih  this  collection.  Twelve 
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resonances  on  U  and  four  resonances  on  U  were  calculated. 

The  number  of  neutrons  obtained  during  fission  through 

. 

resonance  absorption  on  isotope  U  is  c'J.ua  L  xi 


in  which 


v  >  -  <  ?  >. 
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and  summation  is  over  all  the  resonances  of  U  which  are  found 
in  the  group  with  the  number 

When  calculating  reactors  without  reflectors,  syste^  (l)  was 
solved  by  separation  of  the  variables.  In  the  case  of  bizonal 
,  the  method  of  source  integration  was  used  to  find 
the  Eigenvalue  of  the  problem,  while  the  system  of  equations  (l) 
was  solved  by  the  finite-difference  factorization  method  /4/. 
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Table  2 


Group  cross  sections  for  elastic  scattering 
on  hydrogen 
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0.0741 

0.0241 

0,0107 

0.(XX)8 

4 

2.53j2 

0.8299 

0.2694 

0.O875 

0.0390 

0,0027 

0,(4>03 

5 

3.1010 

1.0265 

0.3332 

0.1486 

0,0104 

0.0009 

0.003 

(} 

4.(1001 

1.2987 

0.5791 

0.0403 

0.CO37 

0.0033 

0.0002 

4 

4,7109 

2.1032 

0.1469 

0.0134 

0.0030 

0,0007 

0,0002 

s 

7.2V22 

0.5083 

0.0466 

0.0104 

0.0323 

0.0005 

0.0001 

4.2C40 

0.3903 

0,(»7I 

0.0194 

0,1043 

0,0010 

0,0003 

10 

4,l0t»j 

0.9296 

0.2074 

0.04629 

0.0103 

0.0!  i23 

0,0007 

11 

5,4274 

1,2111 

«;.2702 

0.0603 

0.0134 

0,0030 

00009 

I. 

5,4274 

1.2110 

0.2702 

0.0503 

0.0134 

0.0039 

1.2109 

0.27(G 

0.0303 

0.0173 

14 

5.4272 

1.2110 

0.2702 

0.0776 

U 

5.44(44 

1.2139 

0.3487 

_ 

ICi 

5,5912 

1.6059 

!/ 

7.6771  1 

— 

— 

'  — 

— 

— 
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TABLE  3 


Group  cross  sections  of  non-elastic 
scattering  on^pl^tonlum 


1 _ t 

1 

1  ’ 

3 

1 

4 

5 

6 

2 

0.32 

- 

3 

0.32 

0.32 

— 

4 

0.23 

0.3) 

0.28 

5 

0.09 

0.18 

0.12 

0.15 

6 

0.03 

0.07 

0.04 

0.03 

0.13 

7 

0.01 

0.02 

0.01 

0.01 

0.01 

0.04 

System  of  heat  constants 


TABLE  4 


e  lement 

•A 

Vr 

•ft(l 

PuM 

1026 

746 

2.92 

9.6 

H 

0.33 

12.7 

0 

0 

— 

4.03 

The  multigroup  system  of  constants  for  plutonium,  oxygen  and 

hydrogen  in  given  in  Tables  1-4.  The  corresponding  multigroup 
235  238 

constants  for  U  and  W  are  taken  fo  be  the  same  as  in  the 


paper  by  Marchuk 


mentioned  above.  The  multigroup  water 


of 


the 


constants  used  in  the  calculation  give  a  value 
length  square (^wn  to  indium  resonance,  equal  to  26,.5  cm  ,  3^^=0.665 


•  10  cm  (in  the  diffusion  approximation);  r/\  is  1.'222 


'e.ic 


4 

cm  , 


Results  of  Calculations 
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masses 


The  above  described  method  was  used  to  find  the  critical 

reactors  with  water 


of  spherical 


as  the  moderator.  Both  reactors  without  reflectors  as  well  as 
those  with  a  water  reflector  24..5  cm  thick  were  calculated,  A 
reflector  of  this  kind  is  for  practical  purposes  infinite. 


Figs.  1-7  show'  the  results  of  calculations  for  JHiiMBqpieea^toipS) 


the  active  zone  in  which  is  a  mixture  of  UO  (  •^  =  6  g/cm  )  with 
water.  The  dependence  of  the  critical  masses  of  the  reflectorless 
reactors  on  critical  volumes  at  different  degrees  of  enrichment 
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is  shown  in  Fig.  1,  The  points  with  the  same  ratio  of  the  number 


of  nuclei  *=Ph/Pu«  sire  Joined  by  dotted  lines.  Fig*  2  shows  the 
variation!  in  critical  mass  as  the  concentration  of  the  fuel  in 


the  mixture  changes.  The  correspondih^  dependences  for  reactors 


with  a  reflector  are  given  in  Figs.  3  and  Comparison  of  the* 

Mu 

minifura  critical  masses  of  reactors  with  and  without  a  reflector  are 


given  in  Figure  5.  The  difference  between' 


the 


extrapolated  radius  of  the  reactor  without  a  reflector  and  the 
radius  of  the  active  zone  of  the  corresponding  reactor  wdth  a  re¬ 
flector  is  given  Fig.  6.  For  ref lector less  reactors,  furthermore, 
we  have  calculated  the  cadmium  ratio  in  the  fission  chamber 


CdR 


J  yZ/<f^du. 


Fig.  7  shows  the  CdB.  as  a  function  of  tK  ,  It  corresponds  to  a  1QQ% 

235 

enrichment  v/ith  U 


Similar  calculations  were  also  made  for  the  mixture  UO  -PuO 

2  2 

and  water  (natural  uranium  and  T\jo  ~  ^  g/c®  ) . 


The  results  are  given  in  the  forra  of  graphs  in  Figs.  8-I4. 


The  calculations  were  made  for  different  ratios  ^  —  . 

Fig.  16  shows  a  comparliison  of  the  results  of  the  calculation 

of  critical  masses  and(;^asurements  published  in  /!/.  The  dependence 

••  .  «* 

rf=0.7I 

of  the  extrapolation  length  d  =  0.71  oh'iuel  concentration 

in  the  mixture  is  shown  in  Fig.  16. 

The  value  12,7  given  in  Table  4  is  the  transport  cross  section 
for  single-atom  hydrogen.  This  gave  us  a  certain  exaggeration  of  the 

diffusion  length  in  water,  which  was  offset  by  tw>  low  a  value  for 

2  2 

the  neutron  age  =  26.5  cm  j  -CA  =  30  cm  ) ,  But  in  accordance 

with  the  new  experimental  value  of  the  neutron  age  in  water  (  from 

a  report  given  by  A,  Veynberg  in  the  Institute  of  Atomic  Energy  of 

2 

^  26  cm  ) ,  the  transport 
cross  section  has  to  be  varied  in  accordance  with  the  experimental 


the  Academy  of  Sciences  of  the  USSR,  "U 


value  of  the  diffusion  length.  Nevei^theless,  errrors  in  the  critical 
mass  due  to  this  effect  are  only  slight  over  a  wide  range  of  ratios 


» 
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Fig.  )4.  Critical  masses  of  spherical  reactors  UO^  -  HJ)  with  water  reflector 


as  function  of  - 


10, 


■1 

■ 

i8 

IS 

s 

■*■1 

B 

s 

B 

p 

p 

p 

E 

t 

P 

1 

1 

K 

'S 

a 

H 

i 

■ 

■M 

■I 

■ 

■ 

■ 

1 

■ 

1 

B 

i| 

H 

H 

- 1 

= 

— 

■ 

1 

■ 

■ 

I 

■ 

■ 

B 

1 

1 

1 

■ 

■ 

■ 

■ 

■ 

■ 

■i 

m 

■ 

ii 

— 

■ 

B 

■ 

B 

1 

B 

B 

■ 

M, 

B 

B 

B 

S 

II 

B 

Ii 

■11 

B 

1 

'll 

Wl 

■ 

■1 

■i 

■1 

■1 

■i 

_ 

■ 

■ 

■i 

■1 

■1 

■1 

Bi 

BI 

II 

1 

u  ’•■fy  ‘  ■  ^  ^  »  I  I  »■  1  I  I  I  t.  J 

f  s  to  M  JO  Elftlcffiienf'v  ^ 

»  A 


Fig.  Minimal  critical  masses  of  reactors  UO^  -  HJ): 


1)  v,-ithout  reflector;  2)  7dth  water  reflector 
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Fig.  6.  Economy  of  water  reflector  for  spherical  reactors  UO^  - 


Fig,  7.  Gadmiiim  ratio  as  function  of  «-  Ph/^U**'  • 


0  (O- Ph/Pu*-^)- 
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n  to 
*^3  i 

34  ' 


^  > 

O  * 


to  t  ¥  tito*  g  *  iiJ0\2  k  stto^i  k  to* 
Critical  mass,  1 

Fig.  fi.  Critical  masses  of  spherical  reactors  UO^  - 


PuO£-  H^O 


vfithout  reflectors  at  different  ratios  P  pu/fp«‘ 


The  dashed  lines  join  points  with  the  same  ratio 
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as  a  function 


(/ 
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Critical  volume,  1 


Fig.  10.  Critical  masses  of  spherical  reactors  UO^-PuOg-I^O 
with  a  T/ater  reflector 


Fig.  11.  Critical  masses  of  spherical  reactors  ^0^"  PuO^-  H^O 
Tfith  v^ater  reflector  as  function  or  o» 
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Fig.  12,  Minimal  critical  masses  of  spherical  reactors 


UO^-  PuO^  as  a  function  of  p«  p„;jU 

1)  without  relector j  2)  vrith  water  reflector. 


UO.^— PuOj_- H^O  as  function  of  o-^  /pp,. 


Fig.  Hi.  Cadmium  ratio  as  function  of  •-»h/I*Pu  different  P-Pu/Ppm* 


Critical  volume,  1 


Fig.  15.  Critical  masses  of  spherical  reactors  UO^-  vn.th  a  water  reflector 
Solid  line  shows  theoretical  data  at  lOO-^  enrichment  of  uranium:  o  shows 
experimental  data  at  enrichment  :Tith  uranium. 
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I 


Fig.  16.  Length  of 
spherical  reactors 


extrapolation  d  -  0.71 
UO^  -  H^O  at  a-PH-Pt--- 


for 
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•^’he  critical  parameters  of  the  system  are  determined  hy  the 


Peierls  integral  equation  /!/ 


?(/•)=  \  -^Hr')K(r,r’)dV', 
JUr) 


in  which  lae^ir)  is  the  neutron  distribution  function; 

|:{r)-.  Ar(!,4- 3/4-«c)  is  the  reverse  free  path  length  of  the  neutrons 

(N  Is  the  number  of  nuclei  of  matter  In^^unlt  of  volume,  and 

vy.  Sf  are  the  scattering,  fission  and  absorption  cross  sections); 

Wfi  =  -  "  ^  is  the  reciprocal  mean  neutron  multiplication  factor 

‘  '  >*/  +  », 


per  one  collision  with  a  nucleus  (V  is  the  number  of  neutrons 


uring  one  fission)^. 


For  flsiile  fiagter  <  I ,  for  absorptive  matter  'Y  ^  ^  i  during 
scattering  without  absorption  '^  =  1, 
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The  kernel  ^  determined  by  the  expression 


Kiln  .  ".rM-fi. 


r 


I 

(• 


to 

Integration  is  carried  out  with  respect  ak  the  total  voliime  of 
the  system.  For  a  system  consisting  of  fissile  matter  (medium  1) 
and  an  absorptive  matter  (medium  2),  Eq,  (intakes  the  following  form 


?(0  =  :rr«(r')?(O>^(r.  ?)dr 

T.  J 


-f 


in  which  V  and  V  are  the  volumes  of  the  regions  filled  with 
"1  “2 

fissile  and  absorptive  matter. 

To  obtain  an  approximate  solution  to  Eq.(2^  let  us  single  out 
a  spherical  area  of  volume  V  in  the  center  of  the  system  which  in¬ 
cludes  one  of  the  fissile  media  and  the  aredi^of  the  absorptive  medium 
surroundihg  it.  Let  us  term  the  singled  out  region  the  primary  region, 
and  let  us  call  the  remaining  part  of  the  volume  the  supplementing 
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region,  iet  us  integrate  Eq,^2jwith  respect  to  the  volume  of  the 

X 

primary  region  add  assume  (Jp  (r)  *=  const  .  We  get 


-  i. I dV^  a(0  K {r,  P)  dV‘  +  * (r')K(r,  P)dV' . 

(3) 

in  which  c  =  yi./Yi  • 

Let  us  assume  dV»  =  l^dld^l,  in  which  dll  is  the  solid 

—f 

angle  at  which  the  volume  element  can  be  seen  at  point  r*  from 
point  r.  By  changing  the  order  of  integration,  we  can  transform 
Eq.(3) 


Ti  =  T  J'S*  "*■  + 


(4) 


X  As  shown  in  /2/,  the  error  arising  here  has  a  sign-  favorable 
front  the  point  of  view^of  safety,  and  the  absoibute  value  of  this 
error  decreases  as  the  volume  of  the  primary  region  V  increases. 
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'C 

where  L  is  the  optic  distance  with  respect  to  the  from 

the  point  r  lying  in  the  primary  region  and  the  surface  bounded  by 
the  primary  region  j 

u  =  ^  is  the  optic  distance  with  respect  to  the  ray 

between  the  point  r  and  the  current  point  r*; 
u^  is  the  optic  distance  with  respect  to  the  ray  between 
the  point  r  and  infinity; 

Uj  and  u^  are  the  optic  distances  with  respect**  to 
between  the  point  r  and  a  point  r  lying  in  x  medium  1  or  2, 
respectively,  so  that  the  integrand  is 


when  P  In  med,  , 

I 

I  •  when  r'  in  med,  2. 


Let  US  iMignate  t=  u  -L;  t'»u  -L;  u  -L;  ^  =  u  ^  -1 
it,t* ft^  and  too  are  the  optic  distances  with  respect  to  the  ray 
counted  from  the  surface  of  the  primary  region) , 

Then 

J  («-•  +  cf "*•)  du  =  +  «-*" )  at. 

K*  We  will  use  the  word^  optic* for  the  distances  expressed  in 


lengths  of  neutron  free  path 
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Let  us  replace  t*  and  t"  In  the  Integrand  by  the  mean  values 
of  t*  and  t”  over  the  volume  of  the  primary  region  (the  dependence 
of  these  values  on  the  direction  of  tTl/  is  retained  here).  Let  us 
consider  a  case  inK  which  the  primary  region  is  symmetric  (a  sphere 
surrounded  by  a  layer  of  absorptive) ,  In  this  case,  the  values  of 
the  integrals 


«■••)  du 


do  not  depend  on  the  direction  of  cTL  ,  Eq. (4) assumes  the  form 

I,  =  I  +  “"■>  *  +  T  f  '■‘‘'''X 


tm 
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(5) 


iy  1 

If  IrtM  ^supplementing  region^ 


the  critical  value  yi* 
as 


-  -  AO 

for  tne  primary  regionC  if  c  =  Yi/Ya»*^ _ is  fiKed^-— t  determined 


by  Eq 


\ 

ce~'‘')du. 


U»  ) 


If 


upplmenting  region 


,  but  the  matter  in  the 


primary  region  is  homogeneous  in  composition  (c=l),  the  critical 
value  found  from  Eq.(^5jis 


dV. 


Taking  the  values  found  for  and  "Y 

relationship  for  the  critical  pararaeterf  ^ 
value  c  =  'Y  Iy  f  assumes  the  follov.lng  form 


into  account,  the 
,  given  the  fixed 
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(6) 


From  now  on  we  will  take  ^  and  ^  to  mean  the  Eigenvalues 


for  the  primary  region,  calculated  by  exact  methods, 

Ue 

for  example,  the  improved  diffusion  ]|thod  put  forward  by  Romanov 
in  Ref,  /3/.  f.  <1  ^  4-cr^cXZ^J 

Finally,  to  simplify  Eq.(^6)let  us  use  the  fact  that  when 


avsra« 


igin^t^  and  t"  over  the  volume  V  ^the  middle  part  of  the  cross 
section  of  the  primary  region  greatest  wei^t^/2/. 


Hence  we  v^ill  take  it  approximately  that  t*  and  t”  in  (6)  coincide 
with  and  //(2)  for  a  ray  emerging  from  the  center  of  the 


primary  region. 


Limiting  Cases 


Let  us  consider  the  limiting  eases  which  are  satisfied  by 


Eq.  6; 


1.  Medium  2  is  absent, 


The  relationship  (6)  becomes  a  critical  condition  for  the 
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O 

homogeneous  composition  system  obtained  in  /2/ 

Lzj'  — 

TT- 

(7) 

2.  Media  1  and  2  are  homogeneous  in  condition,  c=l.  Eq.  6 
gives  us  the  critical  relationship  (7)  for  a  homogeneous  composi¬ 
tion  system. 

3.  Homogeneous  system;  the  dimensions  of  bodies  from  fissile 

matter  tend  to  0  when  their  number  is  Increased  to  an  unlimited 

extent  with  retention  of  the  relative  volumes  of  the  first  and 

second  media  (v  and  v  ) .  The  radius  of  the  primary  region  tends 
"1  -2 

to  0,  kHU  hence  t,*-*'©  and  0.  .  The  qualitative  dependence  of 

-> 

the  integrand  in  Eq.  6  on  the  optic  length  t  along  the  Ta.Y'  ili  is 
shoT-Ti  in  Fig.  1. 

At  the  limit,  when  the  increase  in  the  number  of  bodies  and 

the  decrease  in  their  dimensions  and  the  distances  between  them 

are  all  uiuiestBt»kod  (i*e.,  when  the  relative  volumes  v  and  v 

~'l  "2 

are  retained),  the  number  of  discontinuities  in  the  integrand 
function  at  the  intersection  point  of  the  ray  and  the  interfaces 
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of  the  media  per  unit  of  length  will  Increase  to  an  unlimited 
extent.  The  relative  dimensions  of  the  areas  of  continuity  , 
expressed  in  free  path  lengths,  tend  to  k — —  andK — — 

■,W,  -l-  •|V,  -f 


for  segments  of  the  ray  lying  in  <4 


media  1  and  2, 


respectively;^ence  the  integral  of  functions  ^“\t  the 

limit  becomes  the  integral  of  the  continuous  function 


•iVt  +  ’ 


in  v^hich  (X^  and  are  the  reciprocal  neutron  free  path  lengths 
in  media  1  and  2. 


Eq.  6  takes  the  form 


I  -  Tf  = 


in  which 


-  =  - iifS - j._L  _ 

7  Ti  «|V,  +  ^  . 
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coincides  with  the  critical  relationship  (7)  for  a  homogeneous 


system  when^^fche  volume  of  the  primary  region  tends  to  0.  The 
value  of  ^  obtained  is  equal  to  the  reciprocal  mean  nuetron  multi¬ 
plication,  factor  for  a  homogeneous  mixture  of  two  substances  in 


the  proportion V  :v  .  Indeed,  for  a  homogeneous 
'1  '2 


mixture 


It 


in  v/hich 


A 


.  it  is  easy  to  see  that  the  mean  value 


coincides  with  the  value  ^  in  Eq. 

4.  No  supplementing  region.  Eq.  (6)  gives  us  Vi=Vi®.  . 

5.  Medium  2  as  an  Ideally  scattering  substance,  a  =  , 

•m 


f  =  0;  j  e-*'  dt  =  1 


set  a"  ;/ expression  on  the  righijhand 

side  of  6  is  finite.  Since  for  a  medium  in  an  ideally  reflecting 

0  00 

shell  ^  =1  and  =1,  the  expression  in  the  leftpand  side 

of  6  is  finite  at  Y  =1. 

'  i 

6.  Medium  2  as  an  ideally  absorbing  substance,  ^  x.~  • 

,,,  0 

ihe  value  determined  by  the  exact  method  coincides  in  this  case 
with  the  critical  value  'y  for  an  Isolated  body  without  a  shell. 


The  expression  on  the 


rlgh^ 


(since  c  4  0) ,  hence  "Yi'T! 


hand  side  of  the  equation 


114 


vanishes 


7,  The  distance  between  regions  of  fissile  matter  Is  Infinitely 
Increased,  the  density  of  medium  2  remaining  the  same.  Here  ( , 

since  the  radius  of  the  primary  sphere  tends  to  infinity.  In  order 
to  to  keep  Eq.(6)l*  the  left^and  side  finite,  ^  should  imu  tend 

to  When  determined  by  exact  methods,  coincides  with  the 
critical  value  ^  body  in  an  infinite  medium  with  the  constant 

_  It 

Convolution  of  Ea.(^6jfor  System  of  Spheres 

Let  us  apply  Eq.  6  to  a  dystem  of  spheres with  radius  R,  located 

(X 

at  equal  intervals  around  the  volume  of  a  sphere  with  radius  ^  filled 
with  an  absorptive  (or  scattering)  material.  Let  us  single  out 
the  primary  sphere  containing  the  central  sphere  in  the  center  of 
the  system. 

Let  us  divide  up  the  integration  Interval  with  respect  to  the 

beam  (O;  t^o)  into  segments,  in  each  of  which  the  ln<tegrand  function 

in (6)  is  monotonic.  The  boundary  points  of  these  segments  are  the 

points  at  which  the  ray  H-  Intersects  with  the  interfaces  of  media 

1  and  2.  Let  us  designate  the  optical  distances  with  respect  to 

the  ray  between  the  surface  of  the  primary  sphere  and  these  points 

as  t  ,  t*  ,  t  ,  t*  »..»  ti  ,  t^  =  t  ,  in  w'hich  k  is  the  number 
-r  'I  '  1  -k  -te 
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of  sphere5intersected  by  the  ray  (see  Fig,  1^,  Then 


+ce~**')dt=-c  +  (\-e)e'‘‘~-  (l-c)*"'*'  + 
+  (I  -  tf)  r' -  ...+  (1  -  tf)  . 


(9) 

Let  us  Introduce  the  mean  values: 

is  the  mean  optic  distance  of  the  path  in  the  sphere, 

and 

is  the  mean  optic  distance  in  medium  2  with  respect  to 

^  >1 

the  ray  \CL  in  the  surfaces  of  the  two  neighboring  spheres. 

When  the  geometric  progression  in  Eq.  9  has  been  folded  and 
the  mean  lengtlB  have  been  introduced,  Eq»^6^takes  the  form 


Ti  — Ti* 


=  0-^ 


-d-C) 


(I  —  J 


Let  us  supplment  Eq,  10  with  relationships  for  t. 


(10) 

-t/Naiid 


for  a  case  in  which  the  saatterlng  matter  is  uniformly 

distributed  in  the  space  between  the  spheres.  Let  us  represent  the 
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optical  radius  of  the  system  as  t  ^  =  c<^  ,  in  which  (PC  is  the 
mean  probability  of  extracting  neutrons  from  the  beam  per  unit 
of  path  ifcengthllf  there  is  no  medium  2,  then 


in  which  is  the  number  of  spheres  per  unit  of  volume, 

S  is  the  cross  section  of  the  sphere,  and 

^  R)  is  the  transmittance  coefficient  of  the  sphere  for 
a  plane  neutron  beam  determined  by  the  following  relationship  /2/ 


«M)  =  J  exp (- 1) dS  - 


H»  (l-f2»,/?)exp(— 2«i^)  /  j  |\ 

2(«.W 


If  there  are  no  spheres,  the  probability  of  extracting  neutrons 
per  unit  length  may  be  represented  as  the  dioiidaaMe 


where  a  is  the  probability  of  extracting  neutrons  in  medium  2; 

2 

/>S(1— ^(cj/?)]  is  the  probability  of  extracting  neutrons  in  the 
sphere  system  from  medium  2,  identical  in  geometric  size  and  position 
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to  the  spheres  of  fissile  matter. 

The  total  probability  of  extracting  neutrons  from  the  beam 
per  unit  length  is  equal  to 


JT-  aj  +  a,  -  a,  +  /»S[^(«2R)  -  ^(a,/?)l. 


(12) 

The  values  and  (ii+t)  in  Eq.fiojare  determined  by  the 
relationships 


in  which  A»oj/pS  is  the  ^ath  length  expressed  in  the  neutron  free  path 
length  in  medium  2  per  one  intersection  of  the  sphere  by  the  ray  Si, 
^  the  function;  f(x)  is  determined  by  Eq.  (ii)l  We  should  ppint  out 
that  the  function  g(x)  satisfies  the  following  approximate  relation¬ 


ships 
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1»  Qualitative  dependence 
+  f^"'*'on  the  optical 
in  the  direction  of  the  ray  0. 

Medium  ^ (circle 


of  function 
distance  t 

s)  Is  dashed. 
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which  is  equivalent  to 


Finding  the  Critical  j^arameters  Asymptotic 


Diffusion  Method 


Let  us  find  the  critical  relationship  for  a  sphere  from 

aUL 

medium  1  with  radius  R  surrounded  by  a  shew  of  matter  2  with  an 


4^ 


external  radius  r  (isolated  primary  region) .  The  asymptotic 
'"0 

solution  of  the  diffusion  equation  satisfying  the  condition  of 
flnlteness  of  the  neutron  density  is  (0  at  r  =  0,  and  the 
condition  at  which  ?,«  vanishes  on  the  extrapolated  external 
boundary  of  the  shell  r  =  r  +0.71  Y/tc  takes  the  form 

-e  ~o  'i  3. 


— r)  at  7,>1 
•7»‘n«»*(T*)(r,-r)  at  ^^<1^ 
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\t^ 


In  which  the  dependence  determined  by  the  relationships 

=  tan  (k-j-)  at  'Y  ^  1  and  Ic  =  th(k'j^)  at  'y  ^  1/7 A/  (also  see 
Romanov's  paper  in  this  collection). 

Let  us  assume  that  at  the  interface  of  the  two  media  the 
asyaiptotic  solution  ?«!  and  <?.c2  satisfy  the  boundary  condition 
into  accoiint  the  deviation  ?»«  neai  the  boundary  from  the  true 
density  <p(r),  satisfying  Eq.(jl).  In  its  simplest  form  the 
approximate  boundary  condition  derived  by  Davison  /3/  takes  the 
form 


CflSIsii  —  EfiZssi 


at  r=  R. 


Substituting  ?,*  into  the  boundary  relationship,  we  obtain  the 


critical  condition 

—  (1  -M(Ti)^ctg«,*(T,)^l  = 

—  n +«i*<Ti)/?ctha^(T,)(r,-/?)J 
_  *1 

—I*  +  (T»)  ^Ctga,*(Tj(r,  -  ^)1  at 

•* 


Ti  > 


I»<  I- 
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V/e  can  obtain  a  relationship  for  the  critical  parameter  of 


the  sphere  at  a  fixed  ration  c  =  ^ by  replacing  ^  by 
and  by  in  Eq*(jL^.  If  we  assume  ,  Eq.  (13^ 

gives  us  a  relationship  determining  the  critical  parameter 

•f 

for  a  homogeneous  composition  op  the  matter  in  the  ttzzl  primary 
sphere. 

In  practice  |;he  asymptotic  diffusion  method  is  used  for  approxi¬ 
mate  calculations  when  the  radius  of  curvature  of  the  interface  is 
either  large^  or  comparable  with  the  neutron  free  path  length  in  the 
media  (roughly  at  R>0,3c<~^  ).  If  this  condition  is  not  Satisfied, 

the  diffusion  method  may  prove  unreliable. 

x 

Example  1.  Spheres  of  Ou(93.5)  with  radius  R  =  4  cm  are  uni¬ 
formly  spread  through  an  Infinite  space  in  a  scattering  medium  with 
constant  absorption  ^  =  1.1  suid  free  path  lengths  *=  10  cm. 

We  are  to  find  the  critical  distance  between  the  spheres.  The  para- 

-1  XX. 

cm 

isotope 

XX  The  values  of  these  p  irameters  are  found  in  accordance  with 
experimental  data  by  measuring  the  critical  masses  given  in  /5/» 
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meters  for  Ou  (93.5):  are:  0»'74>  =  0.25 

X  Ou  (93.5)  is  metallic  uranium  containing  93.5^ 


Let  us  single  out  the  primary  spherical  region  of  the  greatest 


possible  radius  in  such  a  way  that  it  includes  one  of  the  spheres. 

If  there  are  p  spheres  per  unit  of  volume,  the  radius  of  the 

-i-  5 

primary  sphere  r^  =  d  -  R,  in  which  d  =  p  ^  ±f  the  distance  be- 
tween  the  centers  of  the  spheres.  The  ^roblemTs  such  that  w'e 
should  take  t  =  Cka  in  Eq^^.  The  system  of  transcendental  equa¬ 
tions/6)  and  (13)  can  be  solved  by  the  selection  method  r  .  For 

■"0 

example,  at  r^  =  22.6  cm,  Eq.  (13) gives  us  =  0.553,  ~  0.756, 

and  Eq.(6)ls  satisfied  if  we  assume  that  c,  =  0.769,  instead  of  the 
true  value  c_  =  Jo  "  cm,  ‘^  =  0,535  and  'Y'°  = 

0.683,  Eq.  (6)  is  satisfied  at  e.:=  6.647.  At  in  *  18  cm,  = 

.00  A  ^ 

0,538  and =  0,697,  Eq,^6)is  satisfied  at  c^  =  O.669,  and  so  on. 

By  interpolating  we  find  that  the  true  value  corresp6nds  to 
Tq  =  17,8  cm,  l»e»,  the  critical  distance  between  the  spheres  d  =  r^ 
+  R_  =  81.8  cm.  We  should  point  out  that  in  our  example  the  ratio 
(1  -  —  «-(’j+x)qr)  virtually  coincides  with  unity  in  Eq.  10. 

The  difference  between  this  ratio  and  unity  is  a  small  correction 
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Comparative  table  of  physical  parameters  of  system 


of  spheres  and  homogeneous  fissile  medium 


Fissile  material  medium 


System  of  spheres  of  fissile  material 


Probability  of  fission  per  unit  path  length 


Reciprocal  free  path  length 


•  r  A? (*/+**  +  *<•) 


Probability  of  interaction  with 

(V 

spheres  per  unit  p.th  length 

_ _ _ _ 

Probability  of  extraction  of  neutrons 

from  beam  per  unit  path  length 
/given  uniform  distribution  of  scatt¬ 
ering,  5c  is  dej>ermined  by  (12)/ 


Probability  of  scattering  per  unit  path  length 

-  Nig 


Probability  of  absoiTstion  per  unit  path  length 

«  Ntg 


Breeding  coefficient  per  one  fission 


jS  ^  (y<^f  +  0's) 

breeding  H^oefflolent  for  neutrons  per  one 
collision 

/  _  vcrf  i-  crs 

r  '  ^  cr^ 


(i— 


Coefficient  of  multiplication  by  sphere 
of  neutron  flux  incident  on  sphere  and 


undergoing 


at  least  one 


collision  with  nuclei  of  matter  in  sphere, 

"“(s_ 


/,_JL)/£11;UL+X 

\  T»  /  «  Ti 
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In  -MMaiwof ^oblems  met  with  in  practice.  This  fact  relieves 
of  the  need  to  calculate  l-^^nd  (ti+‘*)cp»  exactly,  when  the  scattering 
matter  Is  nonuniformljf  distrlbute<^n  the  space  between  the  spheres, 
and  when  an  exact  determinations  of  these  values  may  prove  laborious. 
II.  Kinetic  Equation  for  Averaged  Neutron  Density 
The  breeding  of  neutrons  in  a  system  consisting  of  a  large 
number  of  spheres  multiplying  a  stream  of  neutrons  inmpinging 
upon  them  is  similar  to  the  microscopic  process  of  neutron  breeding 
in  a  fissile  medium.  A  parallel  can  be  drawn  betwwen  the  concept 
and  the  physical  value*  in  these  cases  (see  Table), 

On  ecount  of'  the  pMnMM,  the  Boltzmann  kinetic  equation 


can  be  used  approximately  t^escribe  the  neutron  transfer  in  the 
sphere  system,  if  we  replace  the  elementary  constants  by  effective 


microscopic  constants,  in  accordance  with  the  Table,  Here  we  have  to 

- - ^  ‘tuJX>> 

take  Ihto  account  tvpiW /s^erlcal  anomalies  in  our  problem.  The 


first  if  the  nonlsotropic  nature  of  the  angular  distribution  of 
a  neutron  flux  eaMplwirom  the  sphere  when  the  incident,  plane 
neutron  flux  A  similar  phenomenon  iiiim^^6reedAing^eutrons 

in  A  fissile  matter  would  be  the  nonisotropic  nature  of  the  angular 

f 


distribution  of  neutrons 


yHjilAjLA, 


during  fission.  Hence  in  order  to 
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use  the  transfer  equatiohs  for  the  problem  of  a  system  of  bodies 
we  have  to  generalize  the  concept  of  transport  path  for  •  case  of 


nonlsotropic 


neutroi^y^durlng  fission. 


^fN(r,Cl)dCl 

tlfiT  A  is  the  number  of  neutrons  per  unit  of  volume  at 
the  point  r,  whose  directions  of  motion  lie  in  the  element  of  the 
solid  angle  around  0-  ,  the  kinetic  equation  for  \'{r.Q)  takes 

the  following  form 


grad  N  ^,Q)-hvaM(r,  Q)  = 

=  »  Jat  (7,  a-)  ?/(S'  ^  5)  <(Q'  +  , 


(u) 

in  which  ?/(2' 2)  </2' is  the  probable  number  of  secondary  neutrons 
with  a  direction  of  motion  fi,  appearing  during  interaction  between 
the  natter  and  the  i^^^tron  with  a  direction  of  motion  in  the  elementary 
cone  dQ  around  Q,  per  unit  of  its  path  length; 

A  ^(0  is  the  number  of  neutrons  from  the  source  per  unit  of 
volume  and  unit  time. 

If  the  angular  distribution  of  the  neutron^^ter  scattering 
and  fission  is  nonlsotropic,  then 
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^  ^  ^  angular  distribution 

in  which  /,(2' -  2)  and  /,(2' -  2)  are/functions  the  probable 

number  of  neutrons  after  fission  and  scattering,  satisfying  the 

normalization  condition. 


By  integrating  Eq.  14  with  respect  to  dil  ,  we  obtain  the 

continuity  equation 

•divy  (r)  =  (?-«)  vn  (r)  +  q(r). 


(15) 


In  which  «(/■)- is  the  rlj^tr on  density,  and 
jQ)V(r,Q)dQ  — is  the  neutron  flux  density. 

Let  us  multiply  Eq.^4)t)y  il  and  integrate  with  respect  to 
dA  .  Taking  it  into  account  that  the  functions  f^  and  f  ^ 
only  depend  on  the  angle  between  U  and  Q',  ,  w^e  get 


<F  y  2  [2  grad  N  (r,  2))  </2  +  =  0*. 


in  which 


•/r  =  «i  ( 1  —  I*  J  +  «/0  —  I'P/)  +  a,, 

7,=*!  (ad*')/, (2 2') ^2'.  7/=  J(2S')/,(Q2')rfS'; 


if,  and  are  the  mean  cosines  of  angular  distribution 


of  the  neutrons  after  scattering  and  fission. 


The  continui^ 


(l5j  and  (l^  coincide  with  the  relevant 


equations  for  isotropic  neutron  distribution  during  scattering  and 


fission,  provided  the  latter  undergo  substitution  of  c<  by 
*(*  Fi)  +  «/0  >  and  ^  by  ^  ,  In  the  problem  of 


a  system  of  interacting  spheres,  it  is  essentialjk*  to 


Mw/iionj 


onisotropic 


thw  neutron#  ^from  the  spheres  into  account,  since  this 


effect  is  unfavorable  from  the  viewpoint  of  the  safe  handling  of 
fissile  matter  flater  on  we  will  show  that  the  mean  cosine  of  angular 
distribution  when  the  nSttttrons  are  bred  by  a  sphere  is  always  negative). 

The  second  anomaly  in  our  problem,  also  unfavorable  from  the 
viewpoint  of  safety,  is  the  deviation  of  the  true  neutron  density 

— f 

distribution  in  the  scattering  mediuB  from  the  mean  value  n(r) .  If 


the  free  path  length  in  the  rcattering  medium  is  small  compared  with 
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i-yUjLMA 

the  distance  between  the  spheres^  a  neutron  efHtPsimg  from  the 

sphere  has  the  probability  of  going  back  to  the  suie  sphere  after 

scattering  and  breeding  again.  Thus,  some  of  the  neutrons  coming 

dk 

from  the  sphere  cause  an  infinite  train  of  successive  multiplica- 
tions.  Since  the  scattering  of  those  which  had  emerged  and  returned 
to  the  same  sphere  occurs  mainly  in  the  close -lying  layers  of  the 
scattering  medium,  the  phenomenon  leads  to  an  increase  in  the 
neutron  concentration: neat  these  spheres,  compared  with  the  mean 
density  in  the  spaces  between  them  n(r),  contained  in  the  continuity 
oquati'en(l^. 

When  these  two  effects  have  been  taken  into  account,  all  the 
approximate:  methods  of  solving  the  transfer  equationsH(for  example, 
the  asymptotic  diffusion  method)  can  be  applied  to  the  problem  of 
the  system  of  interacting  bodies  ,  Consideration  for  these  effects 
is  made  below.  Let  us  first  look  at  a  case  in  which  the  effects  are 
only  slight. 

Example  2.  Let  us  evaluate  the  critical  number  o:  spheres 

i.  t  f  )  Hi  00 

X  This  relationship  is  often  called  the  generalized  Tick  law  /3/. 
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of  Ou(9.35)of  mass  m  uniformly  distributed  in  the  empty  cavity 

of  a  spherical  reflector  with  rellection  coefficient  (Albedo)  A  ^  0, 

-  0.9.  Volume  of  the  cavity  is  V  .  Let  us  assume  that  there 

"0 

is  no  absorptive  or  scattering  matte^n  the  cavity,  that  the  radius 
of  the  spheresR  is  small^compared  im  with  the  critical  radius  R 
of  an  isolate^ sphere  of  E  Ou  (9.35) • 

Applying  the  asymptotic  diffusion  method,  we  find  the  critical 
relationship.  If  v/e  are  given  the  )!6.bedo  of  the  reflecting  shell. 


the  boundary  condition  on  the  internal  surface  of  the  reflector 
takes  the  forn^/ 

_  gra<»?ae  ^  3  l—A 
•irTae  2  1  + A  ' 


Substituting  the  equation  for  (P  in  the  active  zone  (sphere 

'(as 

system),  we  obtain  the  critical  relationship 

(17) 

in  Vfhlch  a  is  the  radius  of  the  system.  For  the  case  in  point,  in 


which  the 


^bedo  is  close  to  unity,  and  the  coefficient  of  neutron 
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breeding  by  each  sphere  is  considerabljr  less  than  unity^  the 
following  decompositions  hold 


ctga<,Araas 


»l,ka 


•trka 


(18) 


Taking  into  account  that  K 


— a,  we  can  represent  Eq,(^17') 


with  consideration  for 


the  following  form 


X 

Let  us  introduce  the  multiplication  factor  Q  ^which  is  defined 
as  the  ratio^iOf  the  xtx  flux  of  neutrons  leaving  the  sphere  to  the 


flux  of  neutrons  Incident  or  irrespective  of  whether  the  neutrons 

striking  the  wftmm  sphere ^reacted  with  the  matter  or  passed  through 

X  are  7 

it  without  undergoing  any  collisions.  The  values  Q  and  Q  ^related 
as  follows 


(20) 

*- 

As  will  be  shown  1-ter,  Q  can  be  approximately  described  by 
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the  following  interpolation  between  the  extreme  Ik  cases  ^  — ?►  0  and 

R  R 
"0 


(21) 


Using  Eqs.  (2ojand(2^,  and  expressing  all  values  in  Eq.  ^19) 

in  terms  of  the  sphere  mass  m,  the  volume  of  the  system  V  fand 

"O 


the  total  number  otf 
the  critical  condition 


spheres  in  the  system  N  ,  we  get 

“0 


(22) 


Ou. 


in  which  M  is  the  critical  mass  of  the  isolated  sphere  (for  ^(93. ij) 


M  =  51.3  kg  /5/)t  ^  is  the  density  of  the  matter  in  the  m 

,  3 

spheres  (for  Ou(93.5)  is  equal  to  18.8  g/cm  ), 


We  should  point  out  that  the  critical  Eqs.  (19)  and  (2o)  do  not 
include  dependence  on  the  angular  distribution  of  the  neutron  flux 
•MSifiag  from  the  spheres.  Ihe  reason  for  this  is  that  on  account 
of  the  proximity  of  the  albedd  of  thJref lector  A  and  unity,  the 
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i 


mean  neutron  density  varies  slightly  along  the  radius  of  the 

angular 

system,  and  MMMMMMitti  for  the  anlsotroply  of  the/dlstrlbutlon 
of  the  flux  during  breeding  gives  a  correction  for  the  small  value 
of  the  gradient  of  the  mean  neutron  density  (  If  the  latter  Is 
constant,  for  example.  In  an  Inflftlhely-extended  system,  the 
anysotropj^  has  no  effect  at  all  on  the  state  of  the  system). 

Multiplication  of  Instant  Neutron  Stream  bv  Body  Made 
of  Fissile  Matter 

Let  us  derive  relationships  for  the  multiplication  factor 
for  a  plane  beam  of  neutrons  Q  and  mean  cosine ^  of  angular 
distribution  Of  the  flux^^'All  the  elementary  processes 

in  the  body  (scattering  of  neutrons  by  nuclei  and  fission)  are 


conildered  Isotropic,  hence  in  the  kinetic 


&V, 


14^ we  must 


assume /(Q' -  2)  —  -^according  to  normalization  of  the  function  _f. 


tel/ 


We  will  takeptsotroplc  sources  in  Eq.(^14^to  mean  secondary  neutrons 


occurring  through  interaction  of  the  incident  plane  beam  and 

the  matter  .  Integrating  Eq.  (14)  along  the  ray  emerging 

from  point  £  in  the  direction  ^  we  get 

X  The  method  used  here  of  reducing  the  problem  v/ith  anisotropic 
source  to  isotropic  distribution  of  sources  is  described  in  /6/. 
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(23) 


« 

.  I 

4K©iV(r.2,  2,)=Jp(r-/2)«(2o.  r  - /2{exp( -Ja(r- 

.  I 

-  V^ydV\dl  +  j ^(2o.  r - /2) cxpl-j a (7-  rQ)^/'!^/.  (23) 


in  which  is  the  direction  of  motion  of  the  neutrons  in  the 
plane  beamj 

1  is  the  coordinate  along  the  ray  «Ds  . 

We  can  obtain  an  expression  for  the  neutron  flux /(Q,Qi),  emerging 
from  the  body  in  the  direction  by  integrating  Eq. (23) with 
respect  to  the  area  of  the  cross  section  normal  to 


4s/(2,2,)-4av  j“iV(;.2.  Qo)dS=  \^{7)va(?;  2,)X 

V 

t  L 

Xexp(-||arf/'|)rfV  +  |^(r.2,)exp  (-iJarf/'DrfV; 


in  which  ij«fn  is  the  optical  distance  between  point  r^  and  the 

—  ^  ^ 

surface  of  the  body  along  the  ray  A  ,  and  V  is  the  volume 
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of  the  body. 


L 

Let  us  designate  x(Q,  r)^exp(  jjae//'l).  .  If  the  body  Is  impinged 
upon  by  a  single  (per  unit  of  area)  plane  neutron  stream  in  the 
direction  ^0  ,  the  density  of  the  sources  ^(r,Qo)> ,  according  to 
the  definition  given  above,  is  equal  to 


^ (r,  2«)  -  P (r)x(-  2,.  r). 


The  expression  for  the  flux  takes  the  folloring  formy 


For  a  homogeneous  body  of  spherical  shape 

a,)  =  JJ  ^vn  (r.(L»'J2,)J  [r(2'U)J  r’t/rd'J '  4 
+  (2'2j]3c|r.  -(Q'5,)lr*</rrf2', 


(24) 

in  which  (QQo)  is  the  cosine  of  the  angle  between  the  unit  vectors 


^  —9 

and  . 

The  multiplication  coefficient  defines  the  number  of  neutrons 
leaving  the  body  per  one  act  of  interaction  between  the  matter  and 
the  neutrons  striking  the  body  is  equal^  to 


- r  •••  •  -  »  • 


in  which 


^  Jy(Q  (r) X (r)r*dr  +  f  (Q r»dr 

U(ni(r —%)dV  jn(r)r*dr 


(25) 


rie 


ZM- pi'-.  (S'2)|-g-,  -.(r)-p|r, 


The  mean 
the  body  is 


cosine  yU  of  the  neutron  flux  emerging  from 


:  -  (r)  n  (*•)  —  f  Ml^/.  (r)J*r’dr 

Jy(iru,)rfa  Qj‘Ur)r*dr 


(26) 


in  whicH- 


ri(r)=J(‘i'L»)z|r(a'2)J  ,.Ai(r)=J(I>'2)«Ir.(2'2))-^. 


0  »> 

When  deriving  Eq.(26)we  made  use  of  the  following  relationship. 
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valid  for  the  arbltary  function  V(Q'n): 


J  Q<r  CJ'Ll)  j  (2'a)  »F  (ii'S)  dii. 


In  a  limiting  case,  in  the  radius  of  the  sphere  R  is  small 

compared  with  the  critical  radius  R  and  with  the  free  path  length 

"0 

of  the  neutrons,  we  get  v«<srx.  Confining  ourselves  to  the 
first  terms  in  the  expansion  of  equalities  (25)  and  (26)  into  a  power 
series  R,j|  we  get 


Qt  =  *,  I*- 


15 


(27) 


Let  us  consider  another  limiting  Base  in  which  the  radius  of 
the  sphere  R  is  close  to  the  physical  radius  R  .  Let  us  expand 
the  neutron  flux  vn  in  Eq.(2j|jinto  a  series  with  respect  to  the 


^Igen functions  ^*\  of  the  Peierls  homogeneous  (l) . 

Omitting  the  operations  described  in  /6/,  we  get 

vn  (r,  (2'S,)I  -  S  (wl,?/  (r,  (2'2,)J. 

/  ^ 

\va\i  =  J ??/ 1''.  (2'2,)1  X \r,  ~  (2'2,)I  (/P. 
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are 


in  which  xcxfct 


The  Eigen  functions 


f: 


the  Eigenvalues  of  Eq,  (l) (/  =  0>  2,.,.). 

form  an  crH^hrf  and  normal  system 


1=J\ 
i  ^  y* 


Confining  ourselves  in  the  expansion  to  the  two  first  terms 

containing  the  function  (p^  ,  which  is  wt  not  dependent  on  (Q^  , 

*♦  -► 

and  C|?  which  is  proportional  to  the  first  power  (fi'fio)  >  and 
taking  it  into  account  that •,.('■)  =To(^).  »pi (/’)■» 0,we  obtain  the  following 
for  a  sphere  close  to  the  critical  state 

g_  _J _ l^>H<r)x(r)riitry  _ 


4 


(27*) 


Let  us  approximately  replace  the  distribution  of  asymptotic 


density  —  parabolic  dependence  ^ 


in  which 

is  the  extrapolated  critical  radius  of  the  sphere  determined 


Oe 


X 

in  the  asymptotic  diffusion  theory  by  the  relationships 


=  /?,  +  0.71t  =  -j-  ,  tg(^). 


Having  integrated  Eq,(^27*),  we  get 


Or 


feft) 

I  -aiRt' 


(28) 


in  which 


-  -i-  -T- 

♦  All  rfi 
T  rfAnh) 


A, 


I  + 


A, 


^  '  A|  111  +  **)— I’ 

AT,  i  +  (i?,  4.  exp  (-  2/?,); 

Af,  -  I  -  JL  f  J?!L  y  ^  J?a.\« . 

«  I  A,.  ;  •  7  I  Ao,  j  • 

L,  =  2  -{-  (/?!,  4.  A/?*,  4. 4/?,  +  2)  J, 


•■p 


Xexp(-2A,); 


X  Here  and  from  now  on  the  linear  dimensions  are  expressed  in 
neutron  free  path  lengths  C<~*  . 
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transmittance  coefficient  for  a  sphere  of  critical  radius  as 
determined  by  (11). 


of 

3) 


Fig,  2,  Multiplication  coefficient  of  a  sphere  made  of  Ou  as  a  function 

the  relative  sphere  size  R/r^  :  1)  dependence  of  1/Qj  2)  dependence  of  l/(f  ; 
experimental  points  derived  by  Stsiborsld.y  and  Kuvshinov;  R  is  radius  of  the 


sphere;  ^  is  the  critical  radius  of  the  sphere. 
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Interpolating  between  the  two  liirdting  *ases,  we  get  the 


following  approximate  expression  for  Q 


/Co  —  /? 


(29) 


Fig.  2  shows  the  dependence  of  Q(R/Rp)  for  Ov^93»5),  ( ^  “ 
0.74,*  =  0.61)  obtained  by  Eq.(2^.  It  also  gives  the 

dependence  Q  (R/R^)  >  determined  by  Eq.^o),  which  can  be  experimen¬ 
tally  [che c ke d^ i r e c 1 1^.  and  we  plot  the  experimental  points  for 

X 

Q  obtained  by  Stslborskiy  and  Kuvshinov. 

We  should  point  Mmt  out  that  the  simplified  expression (2:^ 

X  frcvy.  ^  J 

vel 


in  example  2  conveys  the  (dependence  Q  (R/R  )  j^qualltatlvely  »■■■■■  ,mj , 
except  that  it  slightly  exa&rates  Q  ,  which  is  an  ' 


from  the  point  of  safety 


AnirsoStrouhv  of  Angular  Distribution 

Let  us  derive  a  relationship  for  the  mean  cosine  of  angular 
distribution  of  neutrons  tmerglng  from  Idw  sphere.  Since  p(po=0  , 
in  the  critical  state,  when  =  1,  the  numerator  in  Eq.  (26) regains 
finite  and  the  medn  cosine  vanishes  on  account  of  the  fact  that  Q 
Thus,  near  the  critical  state  the  anisotropty  disappears. 


But  it  cannot  be  disregarded  completely  since  the  expression  for 


and  contains  the  product  which  remains  finite  in  the 
critical  state.  In  order  to  obtain  an  interpolation  expression 
tor i^GL  ,  let  us  continue  analytically  the  dependence  ot on 
R,  as  determined  by  to  the  supercritical  state  in  vdaich  =  1. 
In  this  limiting  case  becomes  the  follovving  expression 

-Q_ _ L- 


•  ~  *  5  *r(0  r*  rfr  J  Jfi*  (r)  t*dr 


in  v;hlch  P  (r)  is  the  part  of  the  function«piI/’,  (Q',  5o)l.which  does 
—  I  — 

not  depend  on  the  angle.  Solution  of  the  asymptotic  diffusion 
equation  gives  us  the  following  expression  for  ^ 


?I»C 


Let  us  assume  roughly  that 


in  v/hich 

0.717=^; 
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R  is  the  radius  of  the  sphere  in  the  state  with 
"1 


0  =  4.4934  is  the  root  of  equation  B  =  tan  0 


Integration  gives  us 


pQt  = 


et(T) 


in  which 


4  J?,  rf|  ’ 

•  1  ^^1(7)  —  1  ^  ^  . 

i?,  4f|  1(1+*:)-!  ^ 

/f,  =  ^  -  /P,  +  1  -  (R\  +  2/?,  +  1)  exp  (-  2/?,); 

3  .  • 


^•='-f(irr+i(^r= 

/?.,  -  2R\  +  4/?»,  -  5/?»,  +  ^  +  4/?., 

+  lOi?,  +  13/?, •+  y)  exp(-  2/?,). 


» 


Fig,  3.  Dependence Q  on  R/Ro  for  a  spiiere  of 
Oy  (83,5), 

Dashed  line  is  extrapolated  relationship 
to  the  supercritical  state  withX^aa  ij 
o  -  experimental  points  of  B.D.  Stsiborskl 
and  M, I,  Kuvshlnov 


The  Interpolation  expression  satisfying  the  limiting  cases 
for  Eqs,(^27Vn<^  (3^^  takes  the  form 


1  _ 

(— 

-  J1  + 

«  \ 

hQt 

5i(:)/ 

■r 
*  % 


(31) 

Fig,  3  shows  the  dependence  of on  for  Ou(93r5)  (R  j  = 

3.335;  ^  =  0,108)  described  by  Eq.  3l»  and  plots  the  experimental 
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V 


points  found  by  Stsiborskiy  and  Kuvshinov. 


Multiple  Multiplication:!  by  Sphere  of  Neutrons  Reflected  :From 


Scattering  Medium 


Let  us  deternine  the  number  of  neutrons  imltted  by  a  sphere 
per  primary  neutron,  taking  Into  account  the  possibility  of  their 
reflection  from  the  layer  of  scattering  medium  adjoining  the  sphere. 


For  the  sake  of  simplicity  let  us  use  the  concept  of  the  multi¬ 


plication  factor  Q  related  to  Q  by 


(2(j.  Let  us  single 


out  around  the  sphere  xh  a  region  of  the  xte  scattering  medium  for 
the  given  sphere,  the  reflection  coefficient  (yi^lbedo).  of  which  can 
be  equal  to  A.  If  a  neutron  strikes  the  sphere,  Q 

X 

neutrons  leave  it  after  multiplication.  Of  thf se  (l  -  A)Q  leave 

Xr- 

the  bounds  of  the  separated  region,  while  AQ  netJtrons  return  to 

x2 

the  sphere  and,  after  multiplication  produce  AQ  neutrons  of  the 

next  generation,  and  so  on.  Thus,  for  each  primary  neutron  the 

,  (I  -i4)Q* 

Sphere  emits  (1— /l)Q*+(l— 

neutrons.  Calculation  of  the  multiple  multiplication  can  be  reduced 

x  XX 

to  the  replacement  of  Q  by  (1  -  A) Q  /(l  -  AQ  )  in  the  ylqtxtxxi 


critical  equations. 


Jo 


ambiguity  in  selecting  the  dimensions  of 
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reflecting  layer,  we  will  consider  that  A  is  the  yilbedo  of  iim  ik 
scattering  iyxKz  layer  of  infiniteness  thiclcness.  The  error  9l 
involved  in  this  assumption  k*  has  a  sign  which  is  favorable  from 
the  viewpoint  of  safety.  If  the  radius  of  the  sphere  is  large  or 
comparable  with  the  neutron  free  path  length  in  the  reflector,  in 
order  to  find  A  we  can  use  the  asymptotic  diffusion  theory,  in 
which  the  expression  for  the  yiilbedo  for  an  infinite  medium  surround¬ 
ing  a  sphere  takes  the  form  /4/ 


(32  ) 

in  which  R  is  the  radius  of  the  sphere  (in  free  path  lengths  in 
the  reflector) , 

D  =  —  is  the  asymptotic  diffusion  coefficient  in  the 

scattering  iriedlun  (kith  absorption);  the  dependence  of  k  on  the 
absorption  constant  of  the  medium  is  given  by  the  relationship 
k  =  thCk-y)  • 

Let  us  find  an  expression  for  the  /Llbedo  of  the  infinite  re¬ 
flector  in  another  limiting  case,  when  the  radius  of  the  sphere 
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/.u 


is  small  compared  with  the  free  path  length  of  the  neutrons  in  the 
reflector  medium.  Let  us  assume  we  are  given  a  neutron  density- 
distribution  n(r)  in  the  reflector  (r  li  the  distance  from  the 
center  of  the  sphere) .  In  the  element  dV  in  the  neighborhood  of 
point  r,  ~ dV  neutrons  are  scattered  every  second.  Taking  the 
scattering  to  be  isotropic,  we  find  that  (rR)dV  neutrons 


arrive  a±  the  xjdanxfx  ^here  fromfelement^  dV  every  sewond,  in 
y?)  =f exp  1-/01^.  m  - 

which  P.  is  the  distance  between  point  r  and  the 

surface  of  the  splire  in  direction  Q,  Q,and'"^  is  the  solid  angle  at 
which  the  sphere  can  be  seen  from  point  r.  Integrating  with  respect 
to  dV,  we  obtain  an  expression  for  the  flux  lestaurt  on  the  sphere, 


i .  e . , 


J--~jvn(r) (r,  /?) r-dr. 


(33) 


Ihe  neutron  density  ia(r)  can  be  expanded  into  two  components j 


qj\  (r)  is  the  density  of  neutrons  of  direct  origin  wrrivlng  at 


point  r  from  the  sphere  without  colliding  with^scatterer  nuclei, 


and 


2^  M 


the  diffusion  density.  If  the  radius  of  the  sphere  R 
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is  reduced  with  the  neutron  stream  from  sphere  J  remaining  un- 

"0 

-pc 

changed,  the  fluw  of  reflective  neutrons  of  direct  origin,  striking 
the  sphere,  is  reduced  in  proportion  to  R,  while  the  flux  com¬ 
ponent  due  to  scattering  of  the  diffusion  neutrons  is  reduced  in 
proportion  to  R  ,  Hence  at  low  R  it  can  be  taken  that  in  (33^ 

vn(r);ssvn,,(r)  « 


(34) 

Substituting  Eq,  (^34')into  Eq.  obtain  an  expression  for 

the  ^Ibedo 


A 


(35) 

A  comparison  with  the  numerical  calculation  shows  that  Eq^5) 
is  equivalent  with  a  high  degree  of  accuracy  to  the  rela-tlonship 
obtained  by  interpolating  Eq,^35)with  respect  to  the  limiting  cases 
R  ^  0  and  R  — ^  , 
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A7  = 


r  l  +  -|./?Ei(-2/?)  " 

^  _ * 

3  1  +  —^- - 

^  2(1— In2> 


(36  ) 


Equation  36  expresses  the  dependence  of  the  ^I'bedo  of  an 
infinite  medlun  oft  the  radius  of  the  spherical  cavity  R  and  tts 


absorption 


constant  Y  for  a  case  in  which  the  radius  of 


the  cavity  is  small  compared  with  the  length  of  the  free  path,  and 
the  diffusion  approximation  is  inapplicable.  For  intermediate  R_ 
the  ^Ibedo  may  be  found  approximately  by  interpolating  between  the 

3 

dependence  curves  for  A(R)  in  Eqs. (32)and(j6j.  For  =1.1,  this 

interpolation  curve  is  given  in  example  3.  We  should  point  out  that 

'^values 

as  ^  Increases,  the^^iiWK  of  ^  restricting  the  sphere  of  applica¬ 


tion  of  Eqs.  1^3?)  and (36]  are  greatly  increased. 

Example  3..  Spheres  of  Ou(93.5)  are  uniformly  extended  through 

u 

infinite  space  in  a  scattering  medium  with  the  absorption  constant 
=  1.1  and  free  path  length  cx^  =  10  cm.  Let  us  calculate  the 
dependence  of  the  critical  distance  between  spheres  on  the  sphere 


radius  R  by  the  r.ethod  set  forth  in  Section  2.- 


An  infinite  system  is  critical  when  the  averaged  breeding  co- 
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efficient  In  the  system  is  equal  to  unity. 

\*  /  «  1*  ’ 

(37) 

Let  us  express  the  mean  probability  of  extracting  neutrons 
from  the  beam  <X  and  density  of  location  of  spheres  p  in  terms 
of  the  distance  between /spheres  d 


and  let  us  take  into  account  the  correction  for  multiplication  of 
the  neutrons.  The  relationship  for  the  critical  distance  between 
spheres  d  takes  the  following  form 

Th^'  dependence  of  A(^  for  a  medium  with  determined 

by  Eqs.  (32) and  (36)  is  given  in  Fig.  4.  Figww  5  shows  the  dependence 
of  the  critical  distance  between  spheres  d  on  R  (curve  1)  derived 
by  Eq*  (^38j^^  For  comparison  Fig.  5  (curve  2)  shows  the  dependence  d(R) 
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obtained  by  the  method  of  approximate  solution  of  the  integral 


equation  put  forward  in  Section  1 


Pig.  4.  Interpolated  curve  of  the  relationship 
between  the  Albedo  of  the  infinite  medium  and  the 
radius  of  the  spherical  cavity  R  (R  expressed  in 
length  of  free  travel  of  neutrons  in  the  medium, 

#2  *  0*1  cm”l,  y2  ■  1.1). 

1  -  the  curve  calculated  according  to  formula  (36); 

2  -  curve  calculated  according  to  formula  (32). 


Figure  5  shows  that  the  results  obtained  by  the  integral  method 
contain  a  greater  safety  margin  than  by  the  method  given  in  Section 
II.  HVhen  the  radius  of  the  spheres  in  an  infinite  system  is  reduced, 
ixtepex  both  tatepa  independent  curves  d(R)  merge,  since  both 
methods  satisfy  the  limiting  transition  to  a  homogeneous  mixture. 
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Fig.  5.  Relationship 
of  the  critical  distance 
between  spheres  from  Oy  (93{5|)  to  th^ 
radius  for  a  system  of 
spheres  unlformally  dlstrl- ^ 
buted  and  Infinitely 
scattering  medium. 


i  i,c» 


^  ■o  '  • 

('Calculation  for  a  limited  system  of  spheres  is  not  more  com¬ 


plicated  than  in  the  example  considered.  If  the  system  has  a 

spherically  shaped  radius  a,  for  the  critical  relationship  we  only 

....... 

need  breeding  coefficient  (making  a  correction  for 

V 

multiple  multiplication  and  anisotroj^y  of  the  angular  distribution) 

/  /•  V 

equal  to  the  critical  parameter  of  the  system 
which  is  determined  by  exact  oathods  (dependence  of  critical  radius 
on  optic  radius  of  system  given,  for  example,  in  /2/, 

Limits  of  Applicablity  of  Methods 
Both  the  methods  considered  are  applicable,  irrespective  of 
the  relative  dimensions  of  the  bodies,  distances  between  them  and 
free  path  lengths  of  the  neutrons  in  the  media.  Although  attention 
has  been  given  in  the  main  to  a  case  in  which  the  scattering  medium 


t*  uniformly  fills  the  space  betwwen  the  aultiplying  spheres,  both 


methods  permit  generalization  for  more  complex  cases  of  distribution 
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of  the  scattering 


i-satter  (for  example,  scattering  matter 


in  the  form  of  spherical  layers  surrounding  each  sphere  in  the 
system,  or  a  system  consisting  of  spheres  of  two  types,  and  so  on). 
We  should  recall  that  by  a  scattering  medium  we  mean  a  medium  Vi'lth 
^  1, although  all  the  arguments  are  still  walld  if  medium  2  is 
a  fissile  matter  with  ^  d  1,  differing  from  medium  1  in  pro¬ 
perties. 

The  difference  in  the  limits  of  applicability  of  the  two 

I 

methods  is  as  follow.  Method  *•  Is  reliable  when  alculatlng 
systems  consisting  of  a  large  number  of  multiplying  bodies,  w'hen 
the  dimensions  of  the  system  are  large  or  comparable  with  the  mean 
probability  of  extracting  neutrons  from  the  beam  per  unit  length 
during  ±x  interaction  with  the  spheres,  i.e.,  when  we  can  use  the 
kinetic  equation  for  averaged  values.  Method  I  does  not  have  this 
limitation,  since  Eq.(^6)is  valid,  irrespective  of  the  number  of 
multiplying  bodies  in  the  system. 

11 

On  ksxxx  the  other  hand.  Method  t  is  applicable  to  calculation 
of  the  system  of  spheres  in  a  scattering  medium  surrounded  by  an 
external  refl'ctor  (see  example  2),  whereas  Method  1  assumes  that 

TL 

there  are  no  external  reflectors  (in  principle  Method  t  may  be 
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generalized  for  this  case  as  well) 


We  should  point  out  that  when  deriving  critical  relationships 

for  systems  consisting  of  a  large  number  of  multiplying  bodies^ 

it  is  assumed  that  the  bodies  do  not  form  a  regular  geometric 

tta. 

but  on  an  aver  age  p^iformly  distributed  through  the  volume 
of  the  system,  whereas  the  arrangement  of  the  bodies  in  actual 
systems  is  more  often  than  not  geometrically  regular.  The  appll- 
cation  of  these  methods  to  regular  gives  a  greater  safety 

margin  since  the  presence  of  singled-out  directions  in  the  gisM 
means  a  substantial  dtrengthening  of  the  mutual  screening  of  the 
spheres. 
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OF 

EXPERIMENTAL  STUDY  g  INTERACTION  BECTEEN  TUQ  SUBSCRITICAL  REACTORS 

by 

V 

A.  Kamayev,  B,  G.  Dubovskii,  V.  V.  Vavilov,  G.  A.  Popov,  Yu.D, 
Palaraarchuk  and  S,  P.  Ivanov. 


I 


Experiments  with  Homogeneous  Reactors 

Experiments  were  conducted  to  determine  the  critical  state  of 

cylindrical 

a  system  consisting  of  two  *jE±k±«a±j^homogeneous  vtactors  without 
reflectors  in  an  ordinary  air  medium  at  different  distances  between 
them. 


The  two  cylindrical  reactors  in  the  experimental  plant  were 
attached  to  a  frame, one  on  a  movable  base  and  the  other  on  a  carriage 
moving  along  the  frame.  The  distance  between  the  cylinders  could 
be  eiMaprii  from  0  to  120  cm. 

The  walls  of  the  cylinders  were  made  of  rustless  steel  (lkhl8n9t) 

^  s 

1.5  mm  thick.  Each  cylinder  was  fitted  with  an  emergency  rod  which 
dropped  into  the^iddle  of  the  active  7,one  when  the  set  power  level 
was  exceeded.  The  active  zone  was  a  water  solution  of  UO  (NO  ) 


2  3  2 

salt  containing  uranium  with  90%  enrichment. 

It  is  obvious  that  when  the  reactors  interacted,  each  reactor 

was  subwcritlcal,  even  when  the  system  was  critical  as  a  whole. 
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The  experiments  determined  the  connection  "between  the  degree  of 
subcrltical  u  of  each  reactor  (wImi  considered  separately  from 
each  other)  and  the  distance  between  them  at  which  the  reactors 
form  a  critical  system. 

The  experiment  was  carried  out  In  the  following  way.  The 


cylinders  were  placed  right  next  to  each  other  and 

tanoouwly  wl^i  a  water  solution  of  UOP'(NO');  [until  tne  critical 

2  3  21 

state  was  reached.  The  cylinders  were  then  [moved  (Quickly  jkpart. 

The  systems  became  suberltlcal.  In  order  to  achieve  criticality^ 
fresh  amounts  of  solution  were  added  simultaneously  to  both  cylinders. 


The  cylinders  were  separated  and  filled  with  fresh  Klniln  amounts 

of  solution  until  both  cylinders  became  critical,  irrespective  of 

V.i 

each  other,  lie.,  until  a  position  In  which  the  interaction  was 
equal  to  0.  ixxliKKi  The  inspection  and  control  system  was  carried 
out  in  such  a  way  that  it  was  possible  to  follow  the  behavior  of 
each  eactor  separately,  and  the  system  as  a  whole. 


^J^esults  of  exper'ments  with  the  homogeneous  reactors  are 


shown  in  Table  1  and  in  Figs,  1  and  2,  In  each  separate  (^he~ two 


interacting  reactors  were  of  the  same  size  and  had  the  same  active 


3one  composition.  The  effective  £u  1-reeding  factor  of  the  reactor 
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was  calculated  by  the  following  formula 


eff""  +  **s)  (I +»»/.*) 


(1) 


In  which  X  is  the  geometrical  parameter  of  the  a±±±Kx  critical 
"O 

reactor  and, 

X  is  the  geometrical  parameter  of  the  subcritlcal  reactor. 
In  experiments  1,  2  and  3  the  reactors  ^•••-irtydBBwy'ref lec¬ 
tors,  while  in  experiments  A  and  5  the  reactors  were  set  up  on 
a  graphite  base  forming  the  lower  end-reflector  and  a  vertical 
side  wall/ along  which  the  mobile  cylinder  moved.  The  dependence 

of  K  ^  on  the  distance  between  the  reactors,  found  in  the  first 

is 

three  experiments, /shown  in  Fig,  2, 


ESPfirimsnta  yith.Ji?tfirofieagoug  ILeagtor? 


The  experimental  plant  constituted  a  steel  reservoir,  2,5 
meters  high,  2,1  meters  in  diameter,  with  walls  5  mm  thick.  At  the 


bottom  of  the  reservoir  was  a  steel  base  plate  with  fastenings  to 
which  were  attached  two  guiding  grids  made  of 


for  urnaium  block  channels.  The  active  zone  was  assembled  from 


uranium  blocks/^  and  10/^  enrichment.  The  distance  between  the  reactors 
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Results  of  experiments  with  homogeneous  reactors 
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9  is  so  n  m  r.cm 
Fig.  1.  Critical  volume  of  system  of  tVo 
reactors  as  function  of  distance  between 
interacting  reactors:  1)  fourth  experiment 
2)  first  experiment;  3)  fifth  experiment; 
4)  second  experiment;  5)  third  experiment. 


Fig.  2.  K  cc  one  of  homogeneous  interacting 
reactors  al^^function  of  distance  between  them. 
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Fig.  3,  Kg££  of  one  of  the  interacting 

reactors  as  function  of  distance- between 
them;  1)  first  experiment;  2)  second  experiment 
3)  third  experiment. 


Fig.  4.  K  of  one  of  the  interacting  reactors 
as  function  of  distance  between  them: 

1)  first  experiment;  2)  second  experiment;  3)  third  experiment 


was  changed  by  reloading  the  MNMfefa^/ln  the  lattice.^he  critical 
state  was  achieved  by  filling  the  reservoir  with  water.  The 
experimental  method  was  similar  to  that  described  for  homogeneous 
reactors.  Results  of  the  experiments  are  shown  in  Table  2  and  in 
Figs,  3  and  The  experiiiients  demonstrated  that  the  combination 
of  two  roughly  equilateral  subcritical  reactors  flush  with  one 
another  (k  ^  =  0.94)  is  critical  for  each  of  them.  Two  identical 
subcritical  cylindrical  reactors  with  =  1  and  K^^<0,94  each, 

not  exceeding  2,0  at  cannot  form  a  critical  system. 

The  effective  Interaction  of  two  identical  subcritical  reactors 
at  the  given  K  (breeding  coefficient  in  an  infinite  medium)  is 
determined  by  the  solid  angle  between  the  adjoining  surfaces. 

The  effectiveness  of  the  interacticn  of  two  subcritical  reactors 
in  water  at  large  distances  is  considerably  less  than  in  air.  If 
there  is  a  30  cm  protective  layer  of  water  between  the  two  reactors, 
interaction  between  them  can  be  ignored;  here  the  error  in  deter- 
ming  K  will  not  be  more  than  0,1 


fttfcnm.ATTON  OF  SECTIOMLIZED  NUCLEAR  POWER  PLANTS 


by 

G,  I.Marchukpi  B.  G.  Bbbovskly,  V.V.  Smelov  and  Z.  N.  l^ilyntina 


In  connection  with  the  study  of  sectionalized  reactor 

which  a  , 

systems/provide  for  nuclear  energy  in  subcritical  reactors  and 

a  considerable  increase  in  the  burn-op  depth  of  the  fissile 

isotope  /!/,  mathematical  caloolations  of  these  plants  were 

made  and  results  were  obtained  which  can  be  applied  to  certain 

specific  versions  of  the  systems* 

The  physical  system  of  a  sectionalized  reactor  system 

considered  in  this  paper  is  to  some  extent  a  further  development 


and  generalization  of  the  physical  system  of  the  TeSmUk  PUR 
reactor  /2/. 

The  mathematical  calculations  are  based  on  the  use  of 
matrix  spectrorization  of  finite-difference  reactor  equations  /3/. 

1,  We  consider  a  sectionalized  reactor  system  representing 
a  combination  of  a  critical  reactor,  which  is  the  ignition  source 
for \brightening  and  subcritical  sections 

attan^ed  in  series.  On  the  interface  between  the  critical  reactor 


and  the  first  subcrltlcal  section,  and  also  on  the  interfaces 
of  the  other  subcritical  sections  there  are  dirt  layers  which 
only  let  the  neutrons  through  in  one  direction.  This  system  was. 
calculated  in  cylindrical  geometry  for  one  subcritical  section. 

Fig.  1  is  a  sch^atic  of  a  sectionalized  reactor  system. 

In  the  middle  there  is  an  Ignition  critical  reactor,  1,  then 
a  "black"  one  for  thermal  neutrons^  the  barrier, 2,  which  is  a 
combination  of  a  layer  of  uranium  235  of  thickness  0.2  cm,  a 
layer  of  cadmium  and  a  layer  of  [moderatori[>^'ter|  2.5  cm  thick; 
subcritical  section,  3,  adjoins  the  barrier,  and  there  is  a 
"black"  layer  of  uranium  on  its  outside  boundary. 

In  the  uranium  of  the  barrier  layer  tlie  thermal  neutrons 
from  the  critical  Bactor  ■■  are  converted  into  fast  neutrons  of 
the  fission  spectrum,  which  pass  through  the  layers  of  cadmium 
and  water  and  reach  the  eubcritical  section,  where  they  breed. 

The  greater  part  of  the  neutrons  formed  in  the  subcritical  section 
are  unable  to  pass  through  the  neutron  barrier  in  the  reverse 
direction^  since  they  are  moderated  by  the  water  and  absorbed  in 

X 

the  cadmium  layer  . 

x-When  taking  into  account  the  effects  of  the  reverse  transmission 
of  neutrons  by  the  barrier,  as  shown  by  further  calculation,  the 
increase  in  the  effective  breeding  coefficient  of  a  sectionalized 
system  is  not  greater  than  10/s, 
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Fig.  1.  Sectionali^ed  reactor  system:  a)  top  view;  b)  vertlval 
section;  1)  critical  reactor;  2)neutron  barrier; 

3)  first  subscrltical  section. 


2,  The  spatial-power  spectrum  of  the  neutrons  in  each  sub- 
critical  section  which  is  an  element  of  the  sectionalized  subcritical 
system  *an  be  written  in  the  diffusion-age  approximation  as  the 


following  system  of  equations  /3/ 

— =  vOy?  -  S,?  -f  J  ^/r^" 

”  — tT?T  ~  *^l)» 


J-  ■  .  —  '  z=r(u)  ^ _ ^  r 

■  4  2  d/I  ij.  4  2  dn  U.“ 

4-’  2  4  •  2  d/I  15. 


(0 
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Here  *7 -is  the  yield  of  gurtwaxy  secondary  neutrons  per  1 
Liiiillj  Absorbed  neutron, 

x(«)— is  the  energy  spectrum  of  neutrons  occurring  through 
fission, 

S  and  S  are  the  boundaries  of  the  subcritical  section  with 
"I 

precet^ing  and  following  sections,  respectively, 

— ? 

n  is  the  normal  to  the  boundary,  directed  *way  from  the 
center  of  the  sectionalized  system, 

(w)V.  and  t,t  aee  the  flux  densities  /3/  of  fast  and  thermal  neutrons 
(  from  the  preceding  section '-to  the  given  one. 

If  we  use  the  multigroup  representation  /3/,  system  (l)  re¬ 


duces  to  a  system  of  equations  of  the  diffusion  type 
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VDyV»y  -- 
JL  —  JjL 

4  Q  ttn 


---^  V. 


!S. 


?/  ... 

Dj 

dy 

4 

2 

dn  s. 

(y  —  1,  2, . . . ,  m) 


i2) 


If  the  section!  under  consideration  is  close  to  the  critical 
state,  the  iteration  method  of  solving  system (2)  is  not  effective 
In  view  of  the  slow  of  the  solution  /3»4/*  Hence  to 

solve  sys tern (^2)  we  used  method  /3/  based  on  matrix  representation 
of  the  problem  (2) 


—  E‘I»  0, 

4  2  rf«  ii-.  ~  ’  4  2  an  |.v,  *  ’ 

in  which  ^  and  H  are  vectors  with  the  components 
resepctively,  D  and  ^  are  matrices  a  and 


(3) 

(?/)>  {l/l  > 


Sy/ 


-  S/C, 


yt 


-;.OyV^-i-Oy.,^,W7y. 

M.  u  =  lh 
io.  (y  -^  /). 


In  view  of  the  fact  that  the  composition  of  the  barrier  is  a 
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thin:.,  layer  of  uranium  235  and  cadmim,  we  should  discuss  in 
greater  detail  the  formulation  of  diffusion  conditions  on  the 
boundary  between  the  two  media  (A  and  B)  separated  by  a  thin 
absorptive  layer  (C) .  In  the  given  case  it  is  permissible  to 
consider  the  problem  within  the  framework  of  plane  geometry. 

Let  us  take  the  direction  from  medium  A  to  medium  B  as  the 
positive  direction  of  the  normal. 

If  we  use  in  the  plane  problem  to  designate  the 

neutron  flux  in  a  unit  of  volume  of  the  phase  space  (x,Eyit), 
the  integral  neutron  fluxes  y  \x,  |t)  of  each  of  the  energy  groups 
obviously  assume  the  form 

f>+» 


Here  jt^cosB,  els  the  angle  between  the  normal  and  the  direc¬ 
tion  of  motion  of  the  neutron. 

The  diffusion  approximation  corresponds,  as  is  known,  to 

representation  of  the  function  y  (x,  ji)  in  the  form 

*  \ 
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where 


~  -  3£)|i  , 


?yW=  ?/(•*■. 


Assuming  that  the  neutrons  are  not  scattered  in  layer  C, 
but  jcis:  are  only  absorbed  and  reproduced,  we  can  join  the  number 
of  neutrons  entering  layer  C  from  medium  A  every  second  v.'ith  the 
number  of  neutrons  leaving  layer  C  in  the  direction  of  medium  B 
every  second 

'  rf-  ' 

-  3D/;  I"  JiVt*  -- 
0  I» 


rA  I  •! 

K 


y/{*  —  3Da 


ax 


j 


M. 


-r  Q. 
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Here  d  is  the  thickness  of  the  absorbing  layerj 


is  the  macroscopic  absorbtion  cross  section  in  it. 


Q  is  the  nximber  of  neutrons  generated  per  unit  of  time 

X. 

per  cm  of  layer. 


The  balance  of  neutrons  im  flying  in  the  opposite  direction 


can  be  written  in  the  form 


0  0  M 

(*  t*  *  ^  />  4*  *-4 

?£  \  d\t.  —  ZDb  '-r~  )  jiV  —  Q. 

-I  Li 


(5) 

Inte.-ratlng  equalities (4) and (5]  we  get 


-  a 


(6) 


in  which 
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System  (s)  easily  gives  us  the  boundary  conditions  for 
particular  cases: 

a)  Layer  C  does  not  absorb  neutrons  =0): 

”  £ 

In  this  case  the  system  of  equations (6) gives  the  known 
boundary  conditions  on  the  interface  between  the  two  media: 


Oa 


16 


b)  Layer  C  strongly  absorbs  neutrons  (  K  =  )j 

c 

In  this  case  we  get  the  boundary  conditions  on  the  surface 
of  the  "black"  layer: 


I 


•—  ?A  ~  Da 


(6") 


c)  ^edittn  B  is  "black"  for  neutrons  (in  particular,  it  may 
be  i  vacuum)  y. 

Since  it  has  been  assumed  above  that  layer  C  does  not  scatter 

neutrons,  i  neutron  which  has  left  medium  ^  no  longer  returns; 

^ - - 

this  means*  that  layer  ^may|"be  re*arded”^a«a»»ia.  as  a  "black"  one 
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thus  giving  us  the  condition 


6"' 


In  the  matrix  form  wtr  the  conditions  on  the  boundary  between 
two  media  separated  by  a  thin  absorbing  layer  can  be  written  in 
the  form*  of  equations (6j,  in  which  and  are  matrices: 

Ah)  ^ 

while  ^  is  the  matrix  introduced  above. 

The  differential  matrix-vector  equation (3j  aa*  also  all  the 
boundary  conditions  were  represented  in  the  finite-difference 


form,  and  the  solution  of  the  problem  was  found  by  means  of  the 


0^ 


matrix  spectr^lization  method  /3/. 


To  throw  light  on  the  physical  procfesses  occurring  in  a 


sectionalized  reactor  system,  we  carried  out  calculations  of  this 


system. 


V/e  considered  Identical,  homogeneous  lattices  in  the  central 
critical  reactor  and  in  the  annular  critical  section.  The  uranium- 


graphite  heterogeneous  lattice  with  a  20  by  ?0  cm  mesh  consisted 
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of  infinitely  long  cylindrical  rods  of  uraniiim^  3 . 5  cm  in  diameter^ 


with  2%  content  of  uranium  235. 

The  concentration  of  elements  Ihj  the  active  zone  of  the 
3 

lattice  per  1  cm  was  as  follows: 


fi;:-  -  0.(--C2 J".  =  0.0>2061  . 

«■  0.0  H71.  5^  =  0,07077, 

^itjO  ^  0,03346, 


The  thickness  of  the  uranium  235  layer  in  the  neutron 
barrier  ^  -  0,7.  cm,  and  the  thickness  of  the  w.ter  layer  was 
H  =  2.5  cm. 

Calculations  were  made  for  three  values  of  the  effective 


multiplication  coefficient  (K  ,,)  in  the  subcritical  section: 


0.88,  0.93  and  0,97.  In  accordance  with  these  values  of  K  , 


A 

-4 


we  calculated  the  dimensions  of  the  section  AR  =  R_  -  R,, 

(see  Fig.  1).  The  construction  of  the  neutron  barrier  was  the 


sam.e  in  all  cases. 


<4: 


Figs.  2-6  show  graphs  for  the  sp|«ial-energy  distribution 


of  the  neutron  fluxes  in  the  cases  under  consideration. 


Figs.  7  and  8  show  the  fission  integral  Q(r)  for  a  subcritical 
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t. 
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/  - 1(1  •  <i  (f,  W--0.751: 

(r,  w  .=  17.5);  •/  -1')-'  ^  (r. 
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internal  reactor 


The  fission  integral  is  taken-to  mean  the  value 


-  "t 

QW-^  'VI  J  2/rrf«  +  VfT|. 


Fig.  8,  Division  Integral  ^(r)  for  critical 
reactor. 


The  ratio  of  the  maximum  fission  integrals  Q(r)  in  the  suh- 
critical  sections  and  in  the  critical  internal  reactor  for  differaat 
K  ^rr  of  subscritical  sections  is  as  follows 


c 


«,8» 

0,03 

0.07 

.  .  ,  . 

0.73 

0.ft< 

1.21 

/W  A  >  6^  /V 

~  cr 
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As  can  he  seen  from  the  figures  given  in  Figs.  7  and  8, 

/S 

as  in  the  subcritical  reactor  increases  and  approaches  unity, 
there  is  a  considerable  increase  in  the  neutron  flux  in  the  sub- 


critical 


reactor  through  brightening  of  the  high  voltage 


critical  reactor  by  neutrons. 


Thus,  the  calculations  given  confirm  the  initial  assumptions 
that  the  use  of  sectionalized  reactor  systems  makes  it  possible 


to  produce  high  fluxes  of  thermal  neutrons. 
The 


a 


me  ^  _  ”7 

/Batlo  of  ^  to  ^ /where  Q  is  the  mean  Integral  value  of  Jt(r)/ 

max  —  r 

for  subcritical  and  internal  critical  reactors,  and  also  the  ratios 


of  Q  to  Qa  are  as  follows 
^  -c/- 


/C‘«f^0.97 

Critical 

reactor 

.  .  1 ,40 

1.40 

1.40 

2.12 

.  .  .  1,11 

1,.=5«» 

l.W 

1,00  . 
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The  integral  values  of  the  energy  emission  in  the  critical 


reactor  and  subcritical  sections  as  well  as  the  relative  values  of 


the  specific  energy  emissions  per  1  cell  are: 


f-O.O-i 

Critical 

C 

reactor 

-  .  •OMi 

K.  .» 

•K4P 

BIK) 

lOT’ 

Relative  specific 

heat  emlssloa 

per  cell  ....  l.Oa 

i.47 

1.80 

1 

I 


It  follows  from  these  figures  that  the  distribution  of  the 


energy  emission  fna  in  the  subcritical  sections  may  be  Just  as 
advantageous  as  in  the  critical  sections,  and  this  confirms  the 
advisability  of  using  the  sectionalized  reactor  system  proposed 
in  /!/,  It  is  quite  clear  that  if  there  were  a  grpphlte  reactor 
on  the  outer  boundary  of  the  subcritical  sectlonj  instead  of  a 
layer  of  uranium,  dhe  distribution  of  the  endrgy  release  in  th« 


subcritical  section  would  be  still  more  advantageous. 
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